AD- A 143  767 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  ( and  Subtitle) 


Application  of  Surface  Magnetoplasmons 
on  Semiconductor  Substrates 


1 7.  AUTHORf*; 


READ  INSTR 
BEFORE  COMPL 


3.  RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Final  Report 
9/22/80  -  12/31/83 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMB£R('J) 


Donald  M.  Bolle,  A.  V.  Nurmikko,  and  G.S.  Heller  DAAG29-80-K-0074 

jaaPfPUHMTM'U  UHljMNIZ  ATION  NAME  AND  ADDRESS  10.  PROGIUM  ELEMENT.  PROJECT,  TASK 

^  .  AREA  ft  WORK  UNIT  NUMBERS 

--  Bj-nwn  University.  Providence,  R.I.  02912,  and 
Lehigh  University 
Bethlehem,  PA  18015 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 

Post  Office  Box  12211  13.  number  of  pages 

Research  Triangle  Park,  NC  27709 

U.  MONITORING  AGENCY  NAME  ft  ADORESSfi/  drttcrent  from  Controlling  Ottice)  IS.  SECURITY  CLASS,  (ol  this  report) 


12.  REPORT  DATE 

12/31/83 


Unclassified 


1 5e.  DECLASSI  FI  CATION/ DOWN  GRADING 

schedule 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  it  different  from  Report) 


n 


ie.  supplementary  notes 

The  view,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation 


19.  KEY  WOROS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

magnetoplasmons 

substrates 

sub-millimeter  wavelength 
millimeter  wavelength 


20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  identify  by  block  number) 

Attached 


FORM  1 

JAN  7J  14/d  EDITION  OF  1  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 


Abstract 


A  significant  need  exists  ior  developing  active  and  passive- 
components  in  the  millimeter  (KM)  and  subraillimeter  ( EKM )  * 

wavelength,  i.e.,  100-1  GOG  GHz.  The  use  of  solid  state  materials  and 
phenomena,  particularly  ior  semiconductors,  has  become  increasingly 
important  in  this  direction.  Cue  principal  semiconductor  pnenomenon  ® 

whicn  has  been  explored  for  use  in  I-iM  ana  SMH  control  components  is 
the  solid  state  plasma  state.  lerived  from  this  state  are 
controllable  conductive  anc  dielectric  properties  which,  in  turn,  ® 

are  useful  in  controlling  HM  ana  EMI  transmission  ana  propagation. 

Furthermore,  with  the  application  of  a  static  magnetic  field, 

magnetized  solid  state  plasmas  exhibit  a  variety  of  interesting  ana  ^ 

potentially  useful  effects.  Studies  of  electromagnetic  (Em)  wave 

interaction  with  magnetise. a  semiconductor  plasmas  have  revealed  a 

number  of  specific  physical  properties  which  nay  be  used  for  device  * 

purposes. 


In  this  work,  we  study  first  the  propagation  characteristic  of 
EMM  waves  in  one-dimensional  lnsuiute-d  image  guide  structures.  _-The 
dispersion  characteristics  exnioit  non- reciprocity  which  is 
essential  ior  non- reci procai  aevict  uesign.  Feasibility  stucies 
periormea  on  canonical  meueis  pertinent  to  the  design  of  non- 
reciprocal  devices  suen  as  isolators  ana  differential  pnase  shifters 
yieidtd  highly  encouraging  results.  Therefore  we  proceeded  to  an 
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analysis  of  more  practical  structures,  for  example,  the  H-guide 
structures  which  have  two-dimensional  transverse  cross-sections.  It 
was  found  that  only  hybrid  TE-TM  modes  exist  and  true  surface  waves 
propagate  at  frequencies  above  particular  cutoff  frequencies.  The 
effect  of  semiconductor  loss,  dielectric  loading,  and  structural 
parameters  on  the  dispersion  characteristics  have  also  been 
investigated . 

In  addition  to  the  non- reciprocal  waveguiding  property,  GaAs 
materials  aiso  exhibit  significant  optical  nonlinearity  which  allows 
frequency  conversion.  The  nonlinearity  is  enhanced  by  the  presence 
of  electron  plasmas  in  the  medium,  'we  therefore  explored  the  role  of 
solid  state  plasmas  in  nonlinear  mixing  scnemes.  In  particular,  the 
difference  frequency  generation  of  tfl<  waves  by  two  CCf,  lasers  in  a 
multilayer  planar  structure.  The  efficiency  of  power  conversion  is 
approximately  of  the  same  order  of  magnitude  as  that  of  a 
corresponding  semiconductor  bulk  device.  The  effect  of  a  oiasir.g 
magnetic  field  on  the  nonlinear  process  is  also  investigated.  It 
provides  a  means  for  phase  matching  as  well  as  the  enhancement  of 
the  power  conversion  efficiency. 

Finally,  we  discuss  practical  experimental  considerations  which 
apply  to  the  excitation  and  propagation  of  the  surface  plasma  waves 
in  a  basic  one  interface  structure.  Current  materials  status  appears 
to  favor  GaAs  as  the  main  candidate  because  of  the  high  mobilities 
which  have  been  reached  to  date  in  this  semiconductor  at  cryogenic 
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temperatures  (77K).  The  coupling  of  SMM  radiation  to  a  surface 
plasma  waveguide  has  also  been  attempted  by  using  a  far  infrared 
laser;  however,  this  has  been  hampered  by  the  lack  of  spatially 
clean  transverse  laser  beam  structure. 
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1 .  Brief  Review  of  Solid  State  Plasmas 

1.1  Definition 

In  conductors  and  semiconductors  there  exist  large  numbers  of 
free  charge  carriers.  When  the  population  densities  of  such  charge 
carriers  are  sufficiently  large,  a  situation  arises  where  closest 
neighbor  interactions  dominate  the  behavior  of  the  individual  charge 
carriers.  In  addition  to  this  shielding  phenomenon,  the  dense  cloud 
of  charged  particles  also  exhibit  oscillation  phenomena  when  an 
electric  field  is  applied.  The  electrons  and  ions  are  accelerated 
in  opposite  directions  and  tend  to  separate.  However  strong 
electrostatic  restoring  forces  are  set  up  by  this  charge  separation. 
Collective  oscillations  occur,  i.e.,  the  so-called  plasma 
oscillations.  The  word  'plasma’  is  used  to  distinguish  the 
collective  behavior  of  charged  particles  from  the  individual 
behavior  of  the  particles. 

The  shielding  and  oscillation  characteristics  of  plasmas  are 
described  by  two  parameters,  the  'Debye  (or  Thomas-Fermi)  screening 
length  (  ^jj)'  and  the  'plasma  frequency  (  oip)'.  Both  parameters 
depend  crucially  on  the  plasma  density.  In  a  moderately  doped  n- 
GaAs  medium  with,  for  example.  nQ~2. 1 XI 01 ^cm“^ ,  the  plasma  frequency 
is  given  by  (HKS  unit) 

“p2  =  (noe2/me0)  =  10,;Vsec  (1.1.1) 

where 
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ng=  plasma  concentration 

e  =  electron  charge 

* 

m  =  effective  electron  mass 

£q“  permittivity  of  free  space 

The  screening  length  is  given  by 

XD=(2E0e?/3n0e2)1/2  =  230A  (1.1.2) 

where 

Ej.  =  Fermi  energy 

Note  that  for  this  solid  state  plasmas  A  png' ^<1  which  is  very 

different  from  the  condition  Apng'^XXI  imposed  by  the  gaseous 

plasmas.  However,  in  this  range  of  concentration,  the  semiconcuctor 
exhibits  collective  behavior  and  is  still  considered  to  be  in  the 
plasma  state' . 

Another  important  plasma  parameter  is  the  collision  frequency. 
In  an  applied  electric  field,  electrons  or  holes  drift  parallel  to 
the  field.  However,  they  can  lose  their  acquired  drift  velocity 
through  collisions.  If  the  collision  frequency  is  high  compared  to 
the  frequency  of  the  applied  fields,  electrons  do  not  have  enough 
time  to  accelerate  and  decelerate  between  collisions.  Thus,  apart 
from  a  high  frequency  jitter,  there  is  no  separation  of  charges.  The 
system  then  behaves  as  a  lossy  conductor.  On  the  other  hand,  if  the 
collision  frequency  is  small  compared  to  the  frequency  of  the 

applied  fields,  the  charge  separation  must  be  considered  and  the 
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system  is  in  the  plasma  state.  In  doped  semiconducting  media,  the 

collision  frequency  is  determined  by  two  major  scattering 

mechanisms,  namely,  impurity  scattering  and  lattice  phonon 

scattering.  The  probability  of  impurity  scattering  increases  with 

the  concentration  of  ionized  impurities  and  the  probability  of 

phonon  scattering  increases  with  lattice  temperature.  Clearly, 

cooling  will  reduce  the  phonon  scattering  phenomenon.  At  liquid 

nitrogen  temperature  (77°K),  high  quality  GaAs  material  with 
IS  S 

nQ^2.1X10  Jca  J  can  have  a  collision  frequency  in  the  range  of 
10^ 1  Hz.  This  corresponds  to  a  momentum  relaxation  time  of  10-11 
sec  and  an  electron  mobility  of  'ulQ^cm^/V-sec. 

1.2  EM  Vave-Solid  State  Plasma  Interaction 

Consider  an  EM  wave  propagating  inside  a  homogeneous, 
moderately  doped  GaAs  medium.  In  the  MM  and  SMM  wave  range,  the 
interaction  of  conduction  electrons  with  the  applied  EM  field  is 
well  described  by  the  plasma  theory.  The  equation  of  motion  for  a 
free  electron  is  given  by 

qv/ 3t+  (v  .7)  v*  yv  “  (  e/m*  )  ( vX  ( jS  +  B^)  )  ( 1  .2.1  ) 

where 

v_  =  electron  velocity 
V  ■  collision  frequency 

Here  we  have  assumed  that  a  d.c.  magnetic  field  is  applied  to  the 
plasma  medium.  The  above  equation  includes  second  order  terms  such 
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as  the  Lorentz  force  term  v_X_B  and  the  term  (v«v)v.  The  latter  is 
part  of  a  convective  derivative  of  v^  i.e.,  dv/dt  *  5  v/g  t+( v_»V )_v , 
which  gives  the  total  rate  of  change  of  quantities  moving  with  the 
instantaneous  velocity  v^.  These  second  order  terms  (  nonlinear 
effects  )  will  be  discussed  in  more  detail  in  Chapter  3  where  we 
examine  the  use  of  electronic  plasmas  in  nonlinear  generation 
schemes.  Here  we  only  look  for  the  first  order  effects.  Through 
the  equation  that  relates  the  current  density  to  the  carrier 
velocity,  i.e., 


J.  =  -nev  =  o(w)E  (1.2.2) 

the  solution  of  equation  ( 1.2.1)  yields  the  conductivity  tensor. 
This  result  can  also  be  expressed  in  terns  of  a  permittivity  tensor. 
For  a  biasing  magnetic  field  in  the  y-direction,  the  tensor  takes 
the  following  form  : 


0  .2.3) 

(1.2.4) 

(1.2.5) 

(1.2.6) 
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e(0)  •  static  dielectric  constant 
toc  *  cyclotron  frequency 
V  *  collision  frequency 

Note  the  asymmetrical  property  of  this  tensor.  It  characterizes  a 
gyroelectric  material  which,  analogous  to  gyromagnetic  material,  can 
leaa  to  device  designs  which  exhibit  non- reciprocity. 


I 
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2.  Nonreciprocal  Components  Using  Magnetoplasmons  on  Semiconducting 
Substrates 

For  thirty  years,  non- reciprocal  component  such  as  differential 
phase  shifters,  isolators  and  circulators  have  played  a  crucial  role 
in  the  design  of  circuits  and  systems  from  L  band  through  K  band. 
These  several  component  types  have  yielded  systems  with  very 
impressive  performance  characteristics  both  in  discrete  as  well  as 
hybrid  forms  and,  in  recent  years,  been  used  in  integrated 
configurations.  The  most  crucial  aspect  was  the  development  of  very 
high  quality  ferrite  materials  suitable  for  use  at  frequencies  as 
low  as  hundreds  of  MHz  (i.e.,  L  band)  to  as  high  as  the  tens  of  GHz 
in  the  Ku  to  K&  bands.  Rare  eartn  ferrites  and  garnet  materials 
have  allowed  the  design  and  construction  of  devices  with  reasonable, 
even  acceptable,  performance  in  tne  94  GHz  band.  However,  as  we 
move  to  the  SUM  range,  ferrite  technology  will  be  unlikely  to  see 
use  in  the  design  of  non- reciprocal  components  because  of  high  loss. 
Alternative  technologies  are  required  if  we  are  to  have  available 
the  flexibility  and  sophistication  now  possible  using  ferrite  loaded 
devices  at  lower  frequencies.  In  1930,  Nurmikko  and  Eolle-^  proposed 
the  use  of  surface  plasmon  modes  supported  on  appropriately  doped 
high  quality  semiconducting  substrates  as  a  potential  solution. 
These  electromagnetic  modes  coupled  to  the  solid  state  plasma,  and 
subject  to  a  d.c.  magnetizing  field,  provide  a  modal  hierarchy 
exhibiting  strong  non- reciprocal  effects. 


To  design  magnetoplasmon-based  non- reciprocal  devices,  we  need 
to  understand  first  appropriate  dispersion  characteristics.  The 
propagation  of  surface  magnetoplasmons  along  the  interface  between 
two  3emi-infinite  media,  one  of  which  may  be  a  metal  or 
semiconductor,  have  seen  extensive  investigation  both  theoretically 
and  experimentally^-^.  Of  the  many  possible  geometries,  we  shall 
consider  only  the  Voight  configuration  in  which  the  applied  magnetic 
field  is  normal  to  both  the  principal  electric  field  component  and 
the  direction  of  propagation.  Investigations  of  a  class  of  basic 
canonical  structures  have  been  reported.  These  give  insight  into 
the  performance  characteristics  of  elemental  differential  phase 

Q 

shifters  and  non- reciprocal  isolators  .  In  these  early  studies  the 
geometries  considered  consist  of  single  and  multiple  interfaces. 
That  is,  a  single  dielectric-semiconductor  interface  or  slabs  of 
semiconducting  material,  single  or  double  and  loaaed  by  dissimilar 
dielectric  media  on  either  side  and/or  interspersed  between  the 
semiconducting  layers.  Performance  predictions  based  on  a  5-region 
canonical  structure  were  reported'  for  isolators  and  differential 
phase  shifters  in  the  S"M  wave  range.  Sample  results  showed  the 
possibility  of  acceptable  performance  for  isolators  over  a  bandwidth 
of  45  Clin  in  the  5C0  GKs  range  and  of  differential  phase  shifter 
design  over  a  bar.cwicth  of  65  GHs  in  the  }S0  GHz  range.  However, 
for  differential  phase  shifters,  significant  non- reciprocal  phase 
shifts  (= fc!0°-1 20°/ma)  at  sufficiently  low  values  of  insertion  less 


7 


{,<0.76  dB/mm)  can  only  be  obtained  for  improved  relaxation  times 
(  i . e. ,  x=  2.0X1 0-' ^  sec)  than  those  currently  available. 

To  obtain  high  performance  devices  with  available  materials,  we 
explore  new  canonical  structures  as  well  as  optimization  of  design. 
In  section  1  of  this  chapter,  the  propagation  characteristic  of 
magnetoplasmons  in  an  insulated  image  guide  structure  is  studied  and 
the  performance  of  SXM  wave  isolators  and  phase  shifters  is 
evaluated.  To  refine  our  model  and  so  approach  more  closely 
practical  geometrical  structures,  we  examine  in  Section  2  the  H- 
guide  structure  which  has  a  finite  transverse  cross-section.  The 
analysis  and  results  have  also  been  reported  in  references  10  and 
11  . 

2.1  Insulated  Image  Guide  Structures 

2.1.1  Wave  Propagation 

The  structure  to  be  analyzed  is  depicted  in  Fig.  1.  It  is  a 
one-dimensional  structure  which  consists  of  a  magnetized  CaAs  slab 
(region  2)  parallel  to  but  separated  from  a  perfectly  conducting 
plane  by  an  isotropic  dielectric  region  (region  1  ).  In  the  presence 
of  solid  state  plasmas,  the  macroscopic  fields  of  a  propagating  wave 
can  be  obtained  from  Maxwell's  equations: 


vxn  =  -ju)u0H 


(2.1.1) 


VXH  =  j  oj£q£  *E 
V  *  H  =  0 

V*D  =  V*  (eQ£*E)  =  0 

Here,  a  time  dependence  factor  e"  u  has  been  assumed .  Vvo  coupled 
equations  for  the  transverse  electric  and  magnetic  field  £„  ana  1: 
can  be  deduced  from  these  equations,  assuming  a  propagation  factor 
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where 


eg  =  (S2-rf)/S 


(2.1.4) 


All  the  field  components  transverse  to  the  y-axis  can  be  expressed 
in  terms  of  Ey  anf  h'y: 
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where 


y i 2  “  k025+32/3y2 
y22  ’  kQ2ee*32/3y2 
7t  -  x3/3x  -  z 

Since  3/3y*=0,  the  solution  of  eas.  (  2.1.2)  and  (  2.1.3)  can  be 

decomposed  into  pure  TE  or  TM  modes  which  have  nonzero  field 

components  E^,  Hx,  H2  and  Hy,  Ex  and  E, ,  respectively.  Of  the  two 

possible  modes,  only  TM  modes  are  considered  here  since  fields 

corresponding  to  TE  modes  do  not  show  interesting  interactions  with 

the  semiconducting  material.  The  electric  polarization  of  the  TM 

type  of  waves  lies  in  the  plane  transverse  to  the  applied  magnetic 

field.  Therefore,  in  addition  to  the  plasma  oscillation.  the 

electrons  also  process  about  the  applied  magnetic  field.  As  a 

result,  we  get  an  effective  dielectric  constant  c e  which  i3  defined 

by  ea.  (2.1.4).  It  is  a  function  of  the  cyclotron  frequency,  the 

plasma  frequency  and  the  frequency  of  the  propagating  wave.  fig.  2 

shows  a  plot  of  the  effective  dielectric  constant  as  a  function  of 

3  ( 1  )  ( ? ) 

frequency.  The  critical  frequencies  (jJ  ,  Wq  and  define 

regions  where  £g<0.  It  is  in  these  frequency  bands  that  the  bounded 
surface  waves  exist  which  have  exponentially  decaying  fields  in  the 
transverse  direction,  away  from  the  guiding  structure.  We  then,  in 
the  semiconductor  region,  seek  a  solution  of  eq.  (  2.1.3),  together 
with  other  field  components,  in  the  form: 
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where 


and 

k2±  =*  k2  i  jyn/? 
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(2.1.8) 


Similarly, 

in  the  isotropic  dielectric  regions. 
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where 


-k 
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2 


(2.1.12) 


The  unknown  coefficients  A,B,C  and  the  propagation  constant  can  be 
obtained  through  applying  the  appropriate  boundary  conditions,  e.g., 
Hy  and  Ez  must  be  continuous  at  the  x=T  and  x=T+V  planes.  This 
results  in  a  complex  transcendental  equation: 


u  [eek1tanh(k1T)+£1k2  +  ][£  k,  *  £3^'] 

e-2tC2W  »  — - - - — - U - L-±-  (  2.1.13) 

[e  e  *c  1  tanhCk,  T)--  1  k2“][e  ek^  -  c3k2+] 

At  a  given  frequency,  the  complex  propagation  constant  is  obtained 
using  a  HP-1000  computer  and  an  algorithm  based  on  Davidenko’s 
method.  This  algorithm  has  been  reported  in  detail  in  reference 
(2). 


Theoretical  dispersion  curves  have  been  obtained  for  surface 
magneto-plasmons  where  the  effects  of  both  damping  and  retardation 
have  been  included.  Throughout  this  chapter,  calculations  are  shown 
which  use  values  of  parameters  appropriate  for  n-type  CaAs,  i.e., 
e(0)„i3f  lo^/cjp  |=0. 1  ,  v  *0  or  a)  p/60,  and  ujp'101  ^/sec.  Here  we  want 
to  emphasize  the  fact  that  all  frequency  terms  can  be  normalized 
with  respect  to  (jp,  the  plasma  frequency.  All  dielectrics  are 
considered  loss  free  but  the  semiconductor  can  be  considered  either 
loss  free  or  lossy. 

Consider  first  a  single  interface  structure  with  two  semi- 
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infinite  media:  air  and  lossless  GaAs.  The  criterion  for  'true' 
bounded  surface  wave  propagation  along  the  interface  is  that  the 
wave  amplitude  must  decay  exponentially  away  from  the  interface. 
Mathematically,  this  means  that  the  values  of  k2  in  eq.  (  2.1.8)  and 
kj  in  eq.  (  2.1.12)  must  be  real.  These  two  equations  with  k2^=0  and 
k-j  =0  thus  define  dispersion  curves  which  separate  the  surface  wave 
region  from  the  'bulk  magnetoplasmon*  region  and  the  'radiation' 
region,  respectively.  In  the  bulk  magnetoplasma  region,  the 
propagating  waves  become  volume  modal  waves.  The  propagation 
constant  assumes  complex  value  even  if  the  dielectric  loss  is 
neglected  for  both  the  air  and  the  GaAs  regions.  Its  real  part, a  . 
represents  propagation  normal  to  the  surface  which  produces  a  loss 
of  energy,  i.e.,  radiation  from  the  surface  into  the  interior  of  the 
semiconductor.  When  the  semiconductor  loss  is  taken  into 
consideration,  these  volume  modal  waves  suffer  additional 
attenuation  from  its  sinusoidal  field  distribution  within  the  lossy 
GaAs  slab.  In  the  radiation  region,  the  wave  radiates  in  the 
transverse  direction  and  decays  exponentially  in  the  longitudinal 
direction. 

Surface  waves  propagating  in  opposite  directions  experience 
different  dispersion  effects.  This  same  phenomenon  can  also  be 
obtained  by  reversing  the  direction  of  the  applied  magnetic  field. 
This  non- reciprocity  is  manifested  by  the  first  order  linear  term  in 
Y  in  the  transcendental  equation.  This  phenomenon  can  be  utilized 
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for  non- reciprocal  isolator  and  differential  phase  shifter  designs. 
In  Fig.  3  two  branches  were  found  for  both  directions  of 
propagation.  The  lower  branch  approaches  the  light  line 
asymptotically  at  the  low  frequency  end.  At  the  high  frequency  end, 
it  terminates  at  the  intersection  with  the  'bulk*  magnetoplasmon 
dispersion  curve  when  ^ c>0  (  forward  direction  of  propagation  )  and 
it  also  approaches  the  asymptotic  frequency  for  the  unretarded 
surface  magnetoplasmons  when  cjc<0.  This  unretarded  surface 
magnetoplasmon  frequency  is  defined  by 


ua(l,2)-C  w. 


(2.1 .14) 


and  is  sometimes  referred  to  as  the  surface  magnetoplasma  frequency. 
For  the  upper  branch,  whenuic>0,  the  disperoion  curve  begins  at  the 
light  line  at  the  point  where  £  =0  and  it  terminates  when  entering 
into  the  bulk  magnetoplasmon  dispersion  region.  For  u)c<0»  the 
dispersion  curve  begins  at  the  light  line  at  the  point  where r=e^=1 
and  as  the  frequency  increases  it  asymptotically  approach  the 
unretarded  surface  mangetoplasmon  frequency  which,  in  thi3  case,  is 
obtained  from  eq.  (  2.1.14)  by  replacing  £ j  with  e j .  Also  shown  on 
the  same  figure  is  the  surface  magnetoplasmon  dispersion  curve  for 
the  dielectric-GaAs  interface.  There  is  only  one  branch  for  the 
forward  direction  of  propagation,  while  for  the  reverse  direction  of 
propagation  two  branches  are  obtained.  Eackbends  are  observed  when 
damping  is  present1^. 
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When  two  single  interface  structures  are  placed  side  by  side, 
we  obtain  a  single  slab  structure.  The  surface  wave  mode  of  the 
composite  structure  results  from  the  coupling  between  the  individual 

H 

interface  waves  .  For  the  asymmetrically  loaded  single  slab 
structure,  that  is,  a  semiconducting  slab  sided  by  differing 
dielectric  media,  'true'  surface  waves  exist  in  the  w-S  diagram  only 
to  the  right  of  the  light  line  associated  with  the  high  dielectric 
region.  However,  for  some  choice  of  parameters,  we  do  find  the 
existence  of  surface  wave  to  the  left  of  the  high-dielectric  light 
line.  For  example,  in  Figs.  4a  and  4b,  we  find  that  the  S  mode 
exist  between  the  <jJ=1cqC  and  u  “IcqC^  light  lines.  This  is  the  mode 
which  relates  to  the  symmetrical  (  or  even  )  modes  of  a 
symmetrically  loaded  single  slab  structure.  However,  it  only  exists 
for  the  forward  direction  of  propagation.  These  surface  plasmons  are 
called  'generalized'  modes^  since  they  may  have  complex  propagation 
constant  even  when  the  semiconducting  material  is  considered  loss 
free.  Similar  to  the  bulk  magnetoplasmons ,  the  'generalized' 
surface  waves  also  have  a  propagation  component  normal  to  the 
surface  which  carries  energy  away  from  the  interface.  If  we 
increase  the  slab  thickness,  more  energy  radiates  away  from  the 
interface  and  is  trapped  inside  the  semiconducting  medium.  The 
propagating  wave  suffers  from  the  radiation  loss.  When  significant 
material  loss  is  present  in  the  semiconducting  medium,  the  wave 
suffers  sufficient  attenuation  that  it  eventually  vanishes.  On  the 
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other  hand,  if  we  decrease  the  slab  thickness,  more  energy  is  forced 
to  travel  in  the  low  dielectric  region.  The  generalized  surface 
waves  thus  change  into  space  waves  and  we  tend  to  reach  the 
practical  Unit  of  our  numerical  method.  Therefore,  the  generalized 
surface  wave  in  practice  can  be  observed  only  for  a  very  limited 
range  of  the  slab  thickness  and  material  loss. 

The  difficulties  mentioned  above  no  longer  exist  in  the 
insulated  image  guide  structures.  The  presence  of  a  metallic  wall 
on  the  high  dielectric  side  does  allow  a  ' non- radiating'  sinusoidal 
wave  to  exist  in  that  region.  .  Thus  a  true  surface  wave  can  exist 
between  the  light  lines  of  the  low  ana  high  dielectric  regions  on 
the  w-3  plane.  In  Fig.  5a  and  5b  we  display,  as  a  function  of 
frequency,  the  solution  for  the  transcendental  equation  for  both 
directions  of  propagation.  In  these  figures  the  short  broken  lines 
divide  the  surface  magnetoplasma  region  from  the  bulk  magnetoplasma 
region  located  in  the  frequency  ranges  between  f^1 ^  and  fOT  and  above 
f q  .  In  the  surface  magnetoplasma  region  there  exist,  generally, 
three  branches  for  either  direction  of  propagation.  In  the  lower 
frequency  range  we  have,  in  additional  to  the  S  branch,  an  A  branch 
which  relates  to  the  asymmetrical  (  or  odd  )  modes  of  a 
symmetrically  loaded  single  dielectric  slab  structure.  In  the 
higher  frequency  range,  there  exists  a  D  branch  which  describes  an 
electrodynanic  mode  similar  to  the  magnetodynamic  mode  found  on 
ferrite  loaded  wave  guiding  structure.  In  this  frequency  range 
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there  also  may  exist  a  fourth  branch,  the  electrostatic  branch, 
which,  however,  is  heavily  attenuated  and  thus  ignored. 

2.1.2  Non-reciprocal  Isolators 

In  the  design  of  isolators,  the  most  important  parameters  are 
the  isolation  ratio,  the  insertion  loss  and  the  bandwidth.  As  can  be 
seen  in  Fig.  5a  and  5b,  the  A  branch  for  an  insulated  image  guide 
structure  suffers  considerable  attenuation  and  cannot  be  used  in 
practical  applications.  The  S  branch,  however,  exhibits  a  non- 
reciprocal  attenuation  property  which  can  be  used  as  a  basis  for 
isolator  design.  The  isolation  obtained  features  finite  bandwidth. 
On  the  low  frequency  side  it  i3  limited  by  the  decreasing 
attenuation  constants  for  both  directions  of  propagation.  Also, 
little  non-reciprocity  is  observed  here.  On  the  high  frequency  side 
it  is  limited  by  the  increasing  insertion  loss.  This  is  caused  by 
the  fields  gradually  penetrating  into  the  slab  and  becoming 
increasingly  attenuated  at  high  frequencies. 

The  three  primary  parameters  that  affect  device  performance  are 
the  GaAs  slab  width  W,  the  dielectric  slab  thickness  T  and  the 
dielectric  constant  E  ^ .  These  parameters  must  be  chosen  properly  to 
give  the  desired  isolator  characteristics.  However,  increasing  the 
dielectric  thickness  by  a  factor  of  /N  has  approximately  the  same 
effect  on  the  attenuation  and  dispersion  curves  as  increasing  the 
dielectric  constant  by  a  factor  of  N.  Hence,  we  will  only  discuss 


17 


the  effect  on  the  isolator  characteristics  of  varying  the  two 
parameters  W  and  c ^  . 

Figure  6a  illustrates,  for  the  geometry  shown  in  Fig.  1  ,  the 
isolation  ratio  and  insertion  loss  as  functions  of  frequency  for 
various  GaAs  slab  width.  At  a  given  frequency,  a  wider 
semiconducting  slab  increases  the  non-reciprocity  in  the  propagation 
characteristics,  thus  giving  a  higher  isolation  ratio.  At  the  same 
tine,  the  increase  of  the  slab  width  yields  a  higher  insertion  loss 
and  a  trade-off  is  indicated. 

Similarly,  Fig.  6b  shows  that  increasing  the  dielectric  loading 
results  in  principally  increasing  the  insertion  loss,  thus  reducing 
the  useful  bandwidth. 

From  Fig.  6a  and  6b  we  find  that  for  W= 25  ym,  T=25  ym  and  e^=4, 
the  geometry  studied  gives  quite  acceptable  isolator  behavior  with 
an  isolation  ratio  of  10:1  and  an  insertion  loss  of  less  than  0.5 
dB/mm  over  a  bandwidth  of  65  GHz  in  the  400  GHz  range.  This  is  a  15? 
bandwidth  and  is  an  improvement  over  the  result  reported  in  ref. 
(9),  where  the  same  isolation  characteristics  were  obtained  over  a 
10?  bandwidth  for  a  five-region  canonical  structure. 


2.1.3  Differential  Phase  Shifters 

In  Figs.  5a  and  5b  it  shows  that  branch  D,  with  a  proper  choice 
of  the  parameters  ,  is  suitable  for  differential  phase  shifter 
designs  over  wide  bandwidths.  Figure  7  shows  the  differential  phase 
shift  and  insertion  loss  as  functions  of  frequency  for  different 
device  configurations.  The  low  and  high  frequency  limits  are  set  by 
f  and  the  insertion  loss,  respectively. 

OO 

From  Fig.  7  we  see  that  for  £^^  =  13,  T=10  pm  and  W=25  pm,  phase 
shifts  of  35°/mm  to  65°/mm  while  an  insertion  loss  of  less  than  1 
dB/ram  can  be  obtained  over  a  bandwidth  of  ~90  GHz  centered  at  550 
GHz.  The  differential  phase  shift  changes  by  less  than  10;®  over  a  30 
GHz  range,  corresponding  to  a  highly  desirable  slope  of  less  than 
0.5°/GHz. 

In  the  design  of  differential  phase  shifters  we  would  also  like 
to  have  as  large  a  phase  shift/insertion  loss  ratio  as  possible  over 
a  reasonable  bandwidth.  Figure  8  shows  the  ratio  of  differential 
phase  shift  to  insertion  loss  as  a  function  of  frequency  for  several 
device  configurations  of  differing  geometric  dimensions.  Curves  of 
insertion  loss  as  a  function  of  frequency  are  also  included  on  the 
same  graph  to  indicate  the  limits  on  bandwidth  set  by  the  insertion 
loss. 


From  figure  8,  we  obtain  a  differential  phase  shift  per  dB  loss 
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of  about  75°/ <iB  for  an  d,%  bandwidth.  This  is  to  be  compared  with 

the  results  reported  in  ref.  (9)  for  the  five  region  models.  There 

a  ratio  of  1C0°/dB  is  obtained  over  the  same  bandwidth.  However,  a 

- 1  2 

higher  quality  GaAs  material  with  x *20X10  sec  is  required. 

One  final  note  is  that  for  the  two  modes  employed  for  isolation 
and  phase  shifting,  i.e.,  branches  S  and  D,  respectively,  most  of 
the  electromagnetic  energy  i3  guided  along  the  GaAs-air  interface 
and  away  from  the  conducting  plane.  Any  loss  in  the  conductor  would 
have  little  impact  on  the  attenuation  and  dispersion 
cnaracteristics.  The  assumption  of  a  los3  free  conducting  plane  is 
thus  justified. 

2.2  H-guide  Structures 

In  the  previous  section,  we  studied  the  propagation  of  surface 
magneto-plasmon  on  the  insulated  image  guide  structures  and  the 
feasibility  of  using  it  as  a  basis  for  the  non- reciprocal  device 
applications.  The  structures  involved  are  assumed  to  be  uniform  and 
invariant  in  the  y-direction  which  is  parallel  to  the  interface  and 
perpendicular  to  the  direction  of  propagation  z.  The  field 
distributions  of  the  propagating  TM  modal  waves  are  then  uniform  in 
the  y-direction.  In  practice,  however,  the  guided  wave  would  assume 
a  finite  and  non-uniform  transverse  field  distribution,  hence  higher 
order  modes  must  be  considered. 
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Here,  an  investigation  of  the  dispersion  characteristic  is 
presented  for  an  H-guide  structure  (  Fig.  9  )  where  the  two  parallel 
plates  are  taken  to  be  electric  walls  and  the  rectangular  slab  of 
either  infinite  width  (  single-interface  structure  )  or  finite  width 
(  slab  structure  ).  The  slab  is  chosen  to  be  high  quality, 
moderately  doped  n-type  GaAs  material.  It  is  sided  by  two  semi¬ 
infinite  regions  which  have  dielectric  constants  £  j  and  cj, 
respectively.  The  biasing  magnetic  field  is  applied  in  the  y- 
direction  and  is  perpendicular  to  the  electric  walls  and  parallel  to 
the  normal  of  the  sagittal  plane  (  Voight  configuration). 

Electromagnetic  wave  propagation  in  this  structure  is  governed 
by  the  Maxwell's  equations.  In  the  semiconductor  region,  the  field 
distribution  can  be  obtained  from  eos.  (  2.1.2)  to  (  2.1.6). 
However,  for  the  two-dimensional  model,  eq3.  (  2.1.2)  and  (  2.1.?) 
can  no  longer  be  decoupled.  The  presence  of  the  electric  walls 
introduces  coupling  between  the  TE  and  TM  mode  waves  which  greatly 
complicates  the  mathematics  involved.  Fortunately,  the  assumption  of 
perfectly  conducting  planes  enables  the  separation  of  variables.  All 
field  components  have  trigonometric  distributions  along  the  y- 
direction.  For  example,  the  y-compor.ent  of  the  electric  and  magnetic 
field  can  be  expressed  as  follows: 

Ev(x,y)  =  e  (x)  cosk  y 

7  J  J  (2.2.1 ) 

Hy(x,y)  =  hy(  x)  sinkyy 
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where 


ky  *  n^/h  -  transverse  wave  number 
n  ■  integer  =•  order  of  node 
h  ■  height  of  the  structure 

Substitution  of  eq.  (  2.2.1)  into  eqs.  (  2.1.2)  and  (  2.1.3)  yields 
the  following  coupled  linear  equations: 


~  2~  *  ( 2  +Y  2  }  e  i 

3  ^ 


yjjn 


k..h. 


e  y  y 


52b 


3  x 


2  <  Y2^Y2)hy 


•  II— »  l r  0 

5  ye  y 


(2.2.2) 

(2.2.3) 
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2  -  k  2r  -k  2 
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Y  2  ■  lr  2e  _lr  2 

1 2  *0  e  Ky 


The  general  solution  for  the  above  equations  can  be  found  in 

many  texts1^.  The  complex  field  strength  amplitudes  in  the 

semiconducting  region  can  be  written  as  follows: 

For  0$x^ V  Xix  _X1x  X2x  "X2X 

E  =  (e  +Be  +Ce  +De  )cos  k  y 

y  2  y 

X.x  -X  x  X  x  _X-)X 

H  =  q. (e  +Be  )  +  a  (Ce  +De  “  )sin  k  v 

Y  2  1  "2  y 

X  x  -X  x  X?;<  “x-x 

E  .  =  (A  e  +33  e  +CC,e  “  +DD  -a  “  )sin  k  y 

z2  1  1  1  1  y 

XLx  -X  x  x  2X  "X2x 

H  .  =  (An  +E3  o  +CC  a  +DD  e  “  )cos  k  y 

z  2  2  2  2  2  y 

Xlx  “X  ]_:<  X7x  -X2x 

E  .  =  (A  e  +33  e  +CC,o  “  +DD  a  )sin  k  y 

x2  3  3  3  3  '/ 


3  y 

Xlx  “X  x  X  x  -\„x 
=  (A  e  +33,0  +CC  e  2  +DD  a  2  )cos  k  v 

•  •c  +  **  4  4  y  * 
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where 


and 


7  2  l  7  2 

\.2 


KT  =  {(;yi2-^Y22) ‘  +  4kQ2cnky2  )  1/2 


q1>2  =  I  (.CY22-CY12)iRT]/(2toUnk  } 


(2.2.4) 

(2.2.5) 

(2.2.6) 


Here,  { A ^ ,  C^,  D ^ ,  i=1  to  4  are  known  functions  of  y  and  B,  C, 
D  are  unknown  coefficients  of  the  field  distribution. 


If  k^=0,  aa  in  the  one-dimensional  case,  eq.  (  2.2.4)  is 
reduced  to  the  dispersion  equation  of  the  TE  or  TM  modes, 
respectively.  The  normal  modes  of  the  propagating  wave  thus  can  be 
separated  into  pure  TE  or  TM  modes. 


In  the  isotropic  dielectric  regions,  the  field  distribution 
takes  a  simpler  form: 

For  x<0 

E  ,  =  Eek’1'Xcos  k  y 
yl  y 

k-i  x 

H  =  Fe  sin  k  y 
yl  y 
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Again,  F.jJ  ,  i*1  to  A  are  known  functions  of  y  and  E,  F  are 

unknown  coefficients. 


3  For  x>,W,  similar  field  expressions  are  obtained  by  replacing  k1 

1 

,  j  and  £1  with  -k-z  and  £3  and  E1  ,  F1  ,  E,  F  by  another  set  G1  ,  H1  ,  G, 

\ 

j  H. 

* 

!  Matching  the  boundary  conditions  at  x=0  and  V ,  we  can  solve  for 

1  1 

the  complex  propagation  constant  and  the  coefficients  ox'  the  field 
distributions. 

Let  us  consider  first  a  single  interface  H-guide  structure 
consisting  of  semi-infinite  air  and  GaAs  media.  The  lowest  order 
mode  (  k  =0  )  is  TM  in  nature  and  it  has  a  uniform  field 

I  " 

distribution  in  the  y-direction.  Its  dispersion  curve  exhibits  the 

I  "• 

j  same  phenomena  as  those  given  in  Fig.  3  as  well  as  those  retorted  on 

!  5  8 
j  previously^’  . 

|  1 

!  If  the  semiconducting  medium  i3  loss  free,  the  prorogation 

I 

constant  Y  6  for  the  lowest  order  surface  wave  code  i3  purely 


H 


*  4 


»  M 


>  M 


►  m 


1  • 


»  « 


1  « 


1 


24 


> 


« 


imaginary.  It  represents  a  ’true'  bounded  surface  wave  which 
propagates  unattenuated.  The  higher  order  modes,  however,  may  have 
complex  propagation  constant  and  thus  become  generalized  surface 
wave  modes.  Figure  10  shows  a  relationship  between  the  frequency  and 
the  transverse  wave  number  ky.  The  curve  is  defined  by  RT=0  (eq. 
2.2.5  )  which  separates  the  'true'  surface  wave  region  from  the 
'generalized'  surface  wave  region.  For  a  given  frequency,  'true' 
surface  wave  propagation  occurs  for  modes  with  ky  values  outside  the 
shaded  region.  However,  a  very  large  ky/kQ  value  implies  that 
X,  2/k0  and  kj/kQ  also  have  to  be  large  and  the  surface  wave  cannot 
exist  even  within  the  immediate  vicinity  of  the  interface.  Thus  at 
a  given  frequency,  only  those  surface  wave  modes  with  ky  less  than  a 
certain  value  can  propagate  in  the  'generalized'  sense.  And  for  a 
specified  mode,  i.e.,  for  a  given  ky  value,  increasing  attenuation 
occurs  as  the  frequency  decreases  and  eventually  the  mode  becomes 
evanescent. 

The  loss  and  dispersion  curves  for  the  higher  order  surface 
magnetoplasmons  of  a  single-interface  H-guide  structure  are  given  in 
Fig.  11a.  Here  we  have  also  included  the  loss  effect  of  the 
semiconducting  medium.  There  is  only  one  branch  for  either 
direction  of  propagation.  At  the  high  frequency  end,  the  dispersion 
curve  of  the  forward-going  wave  extends  into  the  bulk  magnetoplasma 
region  where  Al f2<0  *hile  the  loss  curve  (long  broken  line)  leads 
into  hign  attenuation  range  where  practical  application  is  limited. 
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For  the  reverse  direction,  the  wave  ceases  to  propagate  as  the 
attenuation  (  indicated  by  the  dotted  line)  becomes  excessive.  This 
is  also  demonstrated  by  the  bendover  of  the  dispersion  curve.  At 
the  low  frequency  end,  note  the  sharp  increase  in  attenuation  below 
cutoff.  The  cutoff  frequency  increases  as  ky  increases.  For  those 
modes  where  the  cutoff  frequency  becomes  equal  to  or  greater  than 
the  surface  magnetoplasma  frequency,  the  surface  wave  ceases  to 

exist.  This  point  is  illustrated  by  a  special  example  shown  in  Fig. 
11b.  In  the  reverse  direction  of  propagation,  the  normalized  cutoff 
frequency  for  the  normal  mode  of  a  single  interface  structure  is 

around  .22  which  is  very  close  to  the  surface  magnetoplasma 
frequency  cj  ^ ^  which  is  shown  in  Fig.  3.  Therefore  wave  propagation 
is  being  cutoff  in  the  reverse  direction  and  the  propagation  becomes 
unidirectional. 

If  the  semiconductor  region  is  thin  enough,  there  will  be 

strong  coupling  between  the  surface  waves  that  propagate  on  both 
surfaces  of  the  slab.  The  H-guide  structure  in  which  we  are 
interested  is  such  a  case.  In  Fig.  11a,  we  use  solid  lines  to 

illustrate  the  dispersion  curve  for  the  symmetrically  loaded,  i.e., 
H-guide  structures.  It  has  two  branches,  S  and  A,  both  of 
which  are  symmetric  about  the  frequency  axis  and  lie  very  close  to 
the  single  interface  dispersion  curve.  As  the  3lab  becomes  even 
thinner,  branch  S  moves  closer  toward  the  light  line,  i.e.,  more 
energy  travels  outside  the  semiconductor  region  and  the  attenuation 


decreases.  However  for  branch  A,  in  the  case  of  a  thinner  slab,  more 
energy  is  trapped  inside  the  semiconducting  slab.  The  dispersion 


curve  moves 

away  from 

the  light  line 

and  the 

attenuation 

increases. 

In  a 

practical 

structure,  the 

semiconductor 

slab 

has 

to  be 

supported 

by  some 

kind  of  substrate. 

Thus 

it 

is 

loaded 

asymmetrically  by  differing  dielectrics  on  either  side.  A  typical 
dispersion  diagram  for  such  a  structure  is  depicted  by  the  solid 
lines  in  Fig.  11b.  Generally  no  bounded  modes  can  exist  in  the 
region  of  theui-3  plane  which  is  to  the  left  of  the  k1  =0  line.  The 
dispersion  curve  has  two  branchs,  S  and  A,  in  the  forward  direction 
of  propagation  and  only  one  branch,  A,  in  the  reverse  direction  of 
propagation.  The  attenuation  for  the  A  branch  is  slightly  lower  for 
those  modes  with  higher  value  of  k  ,  while  for  the  S  branch  it  is 

J 

slightly  higher.  The  attenuation  for  both  modes  can  be  reduced  if 
the  dielectric  constant  of  region  1  is  made  smaller. 

Checking  the  field  distribution,  we  find  a  similar  field 
displacement  effect  for  the  symmetrically  loaded,  single  slab  H- 

Q 

guide  structure  as  that  reported  on  previously  for  its  1- 
dimensional  counter  part.  We  also  find  that  at  low  frequencies,  the 
amplitude  of  the  transverse  electric  field,  Ey,  is  far  greater  than 
that  of  the  longitudinal  field,  E„.  That  is,  the  surface  plasmon  is 
predominantly  of  the  TE  type  and  there  is  little  interaction  with 
the  semiconductor  material.  At  high  frequencies,  however,  Ky>>H„, 
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the  TM  aspect  predominates  and  strong  interaction  with  the 
anisotropic  semiconductor  is  evidenced  by  the  highly  non-reciprocal 
dispersion  curve.  This  non- reciprocity  can  be  used  as  a  basis  for 
the  design  of  isolators  and  differential  phase  shifters. 

In  conclusion,  the  surface  wave  propagation  in  the  insulated 
image  guide  structures  and  H-guide  structures  have  been 
investigated.  Understanding  how  surface  magnetoplasmons  propagate 
along  a  single  dielectric-semiconductor  interface  and  how  two 
surface  magnetoplasmons  on  different  interfaces  couple  to  each 
other,  we  can  readily  predict  the  dispersion  properties  of  more 
complicated  structures  and  obtain  fairly  accurate  assessments  of 
possible  device  performance. 
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5.  Nonlinear  Generation  of  Surface  Plasnons  on  Planar  Structures 

Solid  state  plasmas  exhibit  nonlinear  polarizability  which 
makes  frequency  conversion  possible.  Considerable  interest  exists 
in  the  nonlinear  generation  of  millimeter  and  sutmillimeter  waves  so 
as  to  provide  potentially  tunable  sources  in  the  ICO  to  1 CCO  GHz 
range.  The  nonlinear  generation  method  may  also  be  used  for  the 
excitation  of  surface  magnetoplasmons .  which,  as  was  mentioned  in 
Chapter  2,  can  be  utilized  in  nor.- reci procal  cevice  design.  In  the 
past  twenty  years  tnere  have  been  intensive  efforts  among  scientists 
to  exploit  the  phenomena  of  nonlinear  optics  for  oown  conversion  of 
tiie  signal  frequency  generated  by  tunable  CC2  lasers.  The  CO^  gas 
laser  is  very  appropriate  for  tr.is  application  since  it  combines 
high  power  output  with  many  discrete  frequencies  of  oscillation  in 
the  range  yn.1 1  pm.  Possible  differences  between  two  such 

frequencies  can  range  from  15G-40C0  GHz  while  high  spectral  purity 
with  bandwidth  of  'VICO  KHz  may  be  obtained^’ 


The  basic  theory  of  optical  mixing  in  a  crystal  is  well  known. 
The  interaction  of  propagating  optical  waves  in  a  nonlinear  medium 
induces  a  distribution  of  electric  dipoles  which  radiate  EH  waves  at 
the  sum,  difference  or  nuimonic  frequencies  of  the  input  waves.  In 
tne  dilforence  frequency  generation  case,  this  nonlinear 
polarization  has  a  phase  factor  t-ik„ z)  j_s  determined  by 
the  relative  phase  01  tne  input  waves.  It  is  sometimes  referred  to 
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as  the  'polarization  wave' And  since  this  polarization  wave  will 
generate,  or  drive,  EM  waves  at  the  same  frequency,  it  is  also 
called  a  driving  wave.  In  this  chapter,  the  name  driving  wave  and 
polarization  wave  will  be  used  interchangeably,  as  will  the  name  of 
driven  wave  and  output  wave. 

To  produce  a  favored  output  at  the  difference  frequency  it  is 
necessary  to  have  an  optical  material  with  significant  polarization 
nonlinearity  for  the  input  frequency.  However,  the  nonlinearity  of 
existing  materials  is  not  sufficient  to  give  significant  mixing 
except  over  distances  much  greater  than  the  signal  wavelength.  As  a 
result,  crystals  available  in  relatively  large  sizes  are  preferred 
and  phase  matching,  which  permits  the  coherent  generation  of  output 
wave  over  a  long  interaction  path,  must  be  maintained. 


GaAs  is  one  of  the  most  popular  materials  in  such  an 

application.  It  has  a  relatively  large  nonlinear  coefficient  at 

optical  frequencies  and  small  absorption  loss  at  both  the  optical 

and  SMM  wave  frequencies.  Neglecting  the  effect  of  absorption  and 

boundary  refraction,  the  maximum  output  power  for  the  phase-matching 

1  7 

operation  can  be  obtained  : 
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ui<  2=  frequency  of  the  input  waves 
=  frequency  of  the  output  wave 
2d  =  nonlinear  optical  susceptibility 
L  =  crystal  length  over  which  the  beams  mix 
c  =  light  velocity 

A  =  area  of  the  beam  cross-section 

n1  2  3=refractive  indices  at  2  3 >  respectively 
P,  ■<=input  power  at  the  respective  freauencies 
Rl  2  ^=single  surface  transmission  coefficients 
at  freauencies  to.  9  respectively 
sinc(x,-sin  x/x 

The  efficiency  of  the  power  conversion  is  low.  Co  far,  the  best 
that  has  been  achieved  with  pulsed  input  power  in  the  100  KW  to  1  KW 
range  are  far- infrared  outputs  with  peak  powers  in  the  range  of  1-20 
nVf1y.  CW  operation  was  also  reported  with  two  single  mcae  CW  CC2 
lasers.  With  input  signal  levels  of  25  W  at  10  y  m  wavelengths,  an 

_7 

output  power  level  just  below  10  W  was  obtainea  at  a  frequency  of 
noOCO  GHz  .  Obviously,  the  output  power  at  OHM  wave  would  be  even 
lower  as  a  result  of  the  dependence  of  output  frequency. 

To  raise  the  output  power,  besides  a  higher  input  power,  an 
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increase  of  optical  nonlinearity  and/or  a  decrease  of  the  beam 
cross-section  is  required.  Some  recent  studies  of  harmonic 
generation  at  metal  and  semiconductor  surfaces  and  in  thin  film 
geometries  suggests  that  an  additional  contribution  to  the 
nonlinearity  could  come  from  the  electron  plasmas  present  in  the 
crystal  '  .  Moreover,  the  use  of  thin  film  structures  not  only 

offers  the  desired  reduction  in  beam  cross-section  but  also  leads  in 
the  direction  of  component  integration  at  SMM  wave  frequencies. 

It  is  thus  in  our  interest  to  re-examine  the  process  of 
difference  frequency  generation  in  planar  structures.  The  canonical 
structure  considered  is  given  in  Fig.  12a.  It  has  five  interspersed 
layers  of  dielectric  and  semiconductor  materials.  The  semiconductor 
layers  (the  cross-hatched  regions)  may  be  intrinsic  or  doped  CaAs 
material.  Two  laser  waves,  propagating  collinearly  in  the  z~ 

direction,  are  mainly  confined  within  the  center  layer  which,  in  the 
rest  of  this  chapter,  will  be  referred  to  as  the  optical  waveguide. 
Due  to  the  nonlinearity  in  the  optical  waveguide,  MM  or  SMM  waves 
are  produced  at  the  difference  frequency.  The  propagation 
characteristic  is  determined  by  the  dielectric  property  of  the 
planar  structure  at  that  frequency. 

In  this  work,  special  attention  is  paid  to  the  optical 
nonlinearity,  phase-matching  condition  and  absorption  effect  and 
their  dependence  on  electron  plasma  density  existing  in  the 
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structure.  In  Section  3«1 ,  we  review  the  propagation  characteristic 
of  both  the  optical  and  SMM  waves  in  the  canonical  structure.  In 
Section  3*2,  the  origin  of  the  optical  nonlinearity  in  GaAs  is 
discussed  and  the  distribution  of  the  induced  polarization  source  is 
determined.  In  Section  3*3,  general  expressions  for  the  output 
fields  are  derived  using  two  alternate  approaches.  In  Section  3.4, 
an  explicit  expression  for  the  efficiency  of  power  conversion  is 
obtained  for  a  sample  structure  which  takes  into  consideration  the 
phase-matching  condition  and  the  absorption  effect.  Finally,  in 
Section  3*5,  we  assess  qualitatively  the  effect  of  a  biasing 
magnetic  field  on  the  various  aspects  of  the  difference  frequency 
generation  scheme. 

3.1  Wave  Propagation 

3.1.1  Optical  Wave  Propagation 

Consider  the  propagation  of  optical  waves  in  the  planar 
structure.  In  our  model,  the  optical  waveguide  consists  of  a  thin 
slab  of  GaAs  sided  by  two  thick  layers  of  loss  free  dielectric 
material.  At  optical  frequencies,  the  dielectric  constant,  E^ ,  of 
the  outer  dielectric  regions  should  have  a  much  smaller  value  than 
that  of  the  GaAs  material  so  as  to  ensure  a  good  confinement  of  the 
optical  beams  within  the  GaAs  layer.  On  the  other  hand,  it  should 
have  a  value  large  enough  so  as  to  retard  the  M.M  wave  sufficiently 
to  meet  the  phase-matching  condition.  Quarts,  which  has  a 
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dielectric  constant  of  4  has  proven  to  be  a  good  choice  since  it 
offers  a  compromise  between  the  above  two  requirements.  Also  it  has 
the  additional  advantage  of  a  small  loss  tangent  at  both  the  optical 
and  KM  wave  frequencies. 

For  the  center  GaAs  layer,  its  dispersive  and  absorptive 
characteristics  may  be  changed  by  introducing  free  carriers  into  the 
medium.  However,  at  optical  frequencies,  these  carrier-induced 
dispersion  and  absorption  effects  can  be  neglected  if  the  plasma 
concentration  is  less  than  10lDcn-^.  Then,  we  have  a  dielectric 
waveguide  with  a  thin  slab  of  dielectric  constant  10.8  sandwiched 
between  two  thick  layers  of  lower  dielectric  constant  materials. 
Optical  wave  propagation  in  such  a  structure  is  well  understood  and 
the  analysis  can  be  found  in  many  texts'1^.  Here  we  will  state 
briefly  the  results.  Let  the  wave  propagation  be  along  the  z- 
direction  and  the  structure  uniform  in  the  y-direction.  There  exist 
two  types  of  eigenmodes,  TE  and  TM,  which  have  nonzero  field 
components  Ey,  Hx,  Hz  and  H^,  Ex,  E^  respectively.  Because  of  the 
symmetry  of  the  structure,  the  modes  can  be  further  decomposed  into 
even  and  odd  modes  which  have,  respectively,  a  symmetric  and 
asymmetric  distribution  of  its  transverse  fields  about  the  center 
plane.  The  dispersion  equation  and  field  distribution  of  the  normal 
modes  can  be  obtained  from  the  general  expressions  given  in  Appendix 
1.1  by  setting  tj=t/2  and  t2=u.  A  cutoff  condition  exists  for  the 
nigner  order  modes.  It  is  given  by: 
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2tt t(e-e  1  J1  -  nX  (3.1.2) 

where 

t  *  slab  thickness 

X  *  free  space  wavelength 

e  ■  dielectric  constant  of  the  slab 

■  dielectric  constant  of  the  surrounding  media 
n  *  integer  *  order  of  mode 

With  A~10.6  vim  and  a  slab  thickness  of  2  pm,  all  higher  order  modes 
are  below  cutoff.  Only  the  lowest  order  surface  wave  modes  (TEq, 
TKq)  can  propagate  at  this  frequency.  In  this  work,  the  TEq  mode  is 

preferably  chosen  as  the  input  mode  because  its  fields  are  more 

closely  confined  within  the  slab. 

3.1.2  MM  and  SHM  Wave  Propagatin 

How  we  determine  the  normal  modes  for  the  MM  and  SKM  waves  that 
propagate  in  the  planar  structure.  First  we  have  to  realize  that  in 
this  frequency  range  the  dielectric  property  of  the  extrinsically 
doped  n-type  GaAs  material  depends  strongly  on  the  free  electron 
concentration.  For  our  canonical  structure,  the  outer  GaAs  layers 
are  moderately  doped  with  electron  plasma  concentration  around 
2.1X101 -’cm-^.  The  collision  frequency  is  about  lO^Hz  and  is 
comparable  to  the  MM  wave  frequency.  The  semiconducting  material  is 

in  the  plasma  state.  Its  dielectric  property  is  characterized  by  a 

dielectric  'constant'  defined  by: 

c<„)  -  c(0)-Up2/U2  (3.1.3) 
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It  is  a  function  of  both  the  plasma  frequency  and  the  frequency  of 
the  propagating  wave.  At  MM  wave  frequencies,  the  presence  of 
electron  plasmas  greatly  reduce  the  dielectric  constant.  In  fact, 

1  *■»  _  "7  1  O 

with  Oq-IO  -'em-  ,  the  plasma  frequency  is  in  the  THz,  i.e.,  10  Hz, 
range.  The  doped  GaAs  region  thus  appears  to  have  a  negative 
dielectric  constant. 

The  canonical  model  can  be  simplified  further  by  neglecting  the 
thin  GaAs  layer.  The  simplified  structure  is  depicted  in  Fig.  12b. 
Generally,  this  waveguide  structure  can  support  two  types  of  modes 
TE  and  TM.  The  symmetry  of  the  structure  enables  the  decomposition 
of  the  into  even  and  odd  modes.  The  dispersion  characteristic  and 
field  distribution  of  each  mode  are  again  given  by  the  equations  in 
Appendix  1.1.  In  the  MM  frequency  range  of  interest,  only  the 
lowest  order  TEq  and  TMq  modes  can  propagate  without  being  cutoff. 
However,  there  is  a  major  difference  between  the  propagating  optical 
wave  modes  and  the  MM  wave  modes.  Because  of  the  negative 
-dielectric  constant  involved  at  the  MM  wave  frequencies,  the  fields 
adhere  to  the  interface  instead  of  being  confined  within  the  n-GaAs 
region. 
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3.2  Induced  Nonlinear  Polarization 

It  is  well  known  that  a  dielectric  medium,  when  subject  to  an 
applied  electric  field  E,  tends  to  become  polarized  due  to  the  i 

distortion  of  its  internal  charge  distribution  under  the  influence 
of  the  electric  field.  The  resultant  electric  dipole  moment  per 
unit  volume  is  called  the  electric  polarization  jP.  In  general,  is  • 

related  to  E  through  a  field-dependent  susceptibility  tensor  x(j?) • 

In  i'iKS  units,  this  relationship  is 

P=£0'<(S)E  (3-2.1)  * 

Expanding  X(J-)  as  a  power  series  in  _E,  eq.  (  3.2.0  may  be 

rewritten  as 

» 

P  =  e0[XO0x(2)(E)E]s  (3.2.2) 

where  X  (_E )  is  the  linear  susceptibility  tensor  and  X^^(_E)  is  the  ^ 

quadratic  nonlinear  susceptibility  tensor. 


Similarly,  when  an  EM  wave  is  incident  upon  a  medium,  _P  will 
vary  with  the  _E  component  of  the  wave  in  accordance  with  eq. 
(3.2.2).  The  electric  fields  associated  with  ordinary  monochromatic 
light  sources  are  far  too  small  for  the  observation  of  nonlinear 
phenomena.  Therefore  powerful  optical  lasers  must  be  used.  In  the 
following  two  subsections,  we  will  study  the  contributing  nonlinear 
mechanisms  of  both  the  valence  electrons  and  the  electron  plasmas. 
A  quantitative  estimation  will  be  given  to  the  nonlinear 


» 


» 


» 
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polarization  mechanism  at  the  difference  frequency  by  two  loser 
waves . 

3.2.1  Valence  Electron  Contribution 

CaAs  belongs  to  ?3m  crystal  group.  A  medium  with  such  a 
noncentrosymmetric  crystal  structure  possesses  nonlinear 
susceptibility  of  second  and  higher  orders.  Two  waves  of 

frequencies  and  u) 2  will  then  interact  nonlinearly  in  this  medium. 

In  the  difference  frequency  generation  case,  the  nonlinear 
polarization  at  the  difference  frequency  Au  -  0)2  can  be  written 

in  the  fora: 

=  IB,ye0  (AW  ,U1  ,U2)E6  (-;U1  )ey*(£;w2) 

(3.2.3) 

Here  E  represents  the  total  electric  field  of  the  propagating  wave 
and  a.  B.  Y  refer  to  the  crystal  axes.  The  nonlinear  susceptibility 
y  may  depend  on  all  three  frequencies.  The  symbol  v  is  used  to 

ActBY 

indicate  that  the  contribution  comes  from  the  valence  electrons  and 
the  symbol  *  represents  the  complex  conjugate. 

In  the  KM  and  SKM  wave  range  where  the  interaction  between  the 

EM  fields  and  the  optical  phonons  is  negligible,  the  nonlinear 

1  *7 

susceptibility  is  no  longer  frequency  dependent.  It  becomes  ’(KKS 
units  ) : 


38 


( 


x  =42.9X1Cf12  n/V 


ct*6*y 


otherwise 


(3-2.4) 


Clearly  the  valence  electron  contributions  to  the  nonlinearity  will 
depend  on  the  relative  orientation  of  the  electric  field  with 
respect  to  the  crystal  axes. 


Consider  a  GaAs  waveguide  geometry  in  which  the  surface  normal 
x  coincides  with  the  crystal  axis  a  (001  )  while  the  direction  of 
propagation  z  makes  an  angle  9  with  the  crystal  axis  y(010).  Two 
optical  waves  of  frequencies  and  00  2  propagating  in  such  a 
structure  will  induce  a  nonlinear  polarization  distribution.  At  the 
difference  frequency,  the  nonlinear  polarization  contributed  by  the 
valence  electrons  can  be  expressed  in  terms  of  its  component  in  the 
waveguide  coordinates  (xyz),  i.e.. 


P " 0 ( Au  ) 

vx 

=  e  OX  v  ^  3in29  [Ez(u)  1  )EZ  (u  2  ) -Ey  (w  1  )  Ey  (o>2)] 

+  cos20  [  Ez  (u  1  )  Ey  (w  2  ^  +  Ey  1  ^  Ez  2  ^  ^ 

Pvy^> 

=e0^v^_3in  26[Ex(u1 )E  *(u2)+E  (u, )Ex*(u2)] 


+  c  032  6  [  Ex  ( )  E  z  ( oj  2  )  +  E  z  ( w  1  )  ^2^^ 


pH  WAd) 
vz  '  ' 


,(  sin26[Ex(u)1  )Ez  (  w  2  )  +  E  z  ( u  y  )  z,  ^  (oj2)1 

♦  C0323  [  Ex  (lj  ,  )Ey*(u2)  +  Ey('J1  )Ex*(aJ2)^i 


(3-2.5  ) 
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If  the  crystal  is  cut  in  such  a  way  that  6=45°,  then 

nonlinear  polarization  induced  by  two  TEq  mode  optical  waves 

NL 

only  one  nonzero  component,  i.e.,  the  transverse  component  PyX« 
is  symmetrically  distributed  in  the  transverse  plane  of  the  opt 
waveguide.  Thus  it  gives  a  net,  first  order  contribution  to 
nonlinearly  generated  wave. 

3-2.2  Electron  Plasma  Contribution 

The  presence  of  electron  plasmas  may  increase  the  nonli 
polarizability.  As  was  indicated  in  eq.  (2.2.1),  the  nonli 
effect  comes  in  through  the  Lorentz  force  term  and  the  cor.vec 
derivative  term.  Other  mechanisms  that  result  in  the  second  o 
effect  in  polarization  may  also  exist.  For  example,  the  elec 
plasma  concentration  varies  with  the  pump  field  intensity, 
satisfies  the  following  continuity  equation: 

3n/3t  *  -V-(nv)  (3-2. 

In  the  presence  of  the  pump  fields  of  frequencies  w 1  and  ^ 2 
can  Fourier  expand  all  fields  in  the  form: 


E( r, t) 

*  E(  r;u,  )e"ju1  +  E(  r;;u2)e' jW2t 

♦  E(rvu1^2)e-j(aJrU,2)t  *  E  (  r  2  )  e*  J  (u  1  ""  2  >  t 


♦  c  .  c  . 
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B( r. t) 

-  _B(  r_;w  1  )  e"  1  t  ♦  J3  ( _r;  w  2  )  e  “  2  t 

♦  B(r;^1 +tu2)e- *  B  (  r  ;u>  1  -a.  2  )  e~  *  (w  1  _u  2  }  t 

+  .  +  c.c. 

v(_r,  t ) 

'  v_(L’U  1  )  e”  ^  1  t  +  v(r;ta2)e‘',(J2t 

♦  v(r;;i)1>u)2)e-j(a)1+a)2)t  -  v  (  r  ;w  1  -w  2  )  e~  ^  1  2  > t 

+  .  +  c.c. 

n(_r ,  t ) 

■  ng  +  n(_r;w  ^  )e“  ^  +  n(_r;u)2  )  e_ -^2  *• 

+  n(_r;u  1 2  )  e”  1 +tJ2  ^  t  +  n(_r;u)  i  -oj2  )  e  1  ^2^t 

♦  .  ♦  c.c. 

(3.2.7) 

where  c.c.  represents  complex  conjugate  terns. 

The  only  d.c.  tern  included  in  these  expressions  is  the  nQ  tern 
which  represents  the  steady  state  electron  concentration.  All  the 
other  d.c.  terns  are  omitted  since  these  give  only  higher  order 
contributions  to  the  polarization  nonlinearity. 

Under  the  snail  perturbation  assumption,  we  can  solve  eqs. 
(1.2.1 )  and  (  3.2.6)  by  successive  approximations.  Substituting  eq. 
(  3.2.7)  into  the  above  equations  and  collecting  the  first  order 
terns,  we  get 

v^'^(co)  ■  (je/n  u)E(u)  )  (3.2.3) 
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Throughout  the  rest  of  the  chapter,  the  spatial  dependence  of  all 
the  field  related  quantities  is  implied. 


Collecting  the  second  order  term  of  difference  frequency  An)  ,  we 


get 


v(2)(  Aw) 
j  e/m* 

I - [BUw)*vU  )(  )XB*(U,)] 

A  u)  e- 


1  ^  (oj1  )  *v]_v^  1  ^  ( cup  ) 


(Aui  ) ' 


I  ♦  (1  2  )* 


(3.2.9) 


In  the  above  expression,  (l  H  2)  means  that  similar  terms  can  be 
obtained  by  the  interchange  of  the  roles  of  1  and  2  followed  by  a 
conjugate  operation. 


Now  we  can  derive  the  expression  for  the  nonlinear  current.  At 
the  difference  frequency,  the  current  induced  by  electron  plasmas  is 
given  by: 

«j(Aw  )  *  -e(  Dgv^2  ^  (Acs  )  +  n(uj  1  )_v  ^  1  ^  (in  2  )  +  n  (w  2  )  v/  1  ^  (w  \  )  ) 

(3.2.10) 

which  consists  of  a  linear  part  and  a  nonlinear  part.  If  the 
frequencies  of  the  input  waves  are  much  greater  than  the  difference 
frequency,  i.e.,  when  oj  1  ,  oo 2 ^ >^.to  >  we  get  a  simplified  expression  for 
the  nonlinear  current: 
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» 


JNL(Aw  ) 

to 


JeOwp2e//m*  j  v[E(cj  1  )  •  E*(u  2  )  ] 


to  to 
1  2 


A  to 


(102)]  +to  2  (to  2  )  [  V  •  E.(oo  1  )  1 


to  ^  to  2 


(3.2.11 ) 


In  this  equation  the  last  term  on  the  right  hand  side  is 
significant  only  at  the  boundary  where  the  gradient  of  the  electric 
field  is  large.  In  the  present  analysis,  this  surface  current  term 
is  neglected  and  only  the  bulk  current  is  taken  into  consideration. 


» 


» 


Instead  of  using  the  bulk  current,  the  nonlinear  effect  can 
also  be  expressed  in  terns  of  a  nonlinear  polarisation.  The  bulk 
current  expressed  by  eq.  (3.2.11)  is  equivalent  to  a  polarisation 
term: 


pNL 
— c 


(A<o  ) 


O^p 


2e/m* 


“>  1to2(Ato  )2 


■V  [  JS  (to  j  )  •  E.  (to  2  ^  -1 


(3.2.12) 


where  c  refers  to  the  conduction  electron  contribution.  This  is  a 

general  expression  for  the  nonlinear  polarization  induced  in  the 

interior  of  a  nonlinear  medium  by  two  pump  fields.  With  two  TEq 

modal  waves  as  input,  the  nonlinear  polarisation  induced  through  the 

presence  of  electron  plasmas  has  two  nonzero  components,  the 

transverse  component  and  the  longitudinal  component  P‘^,  both  of 

cx  cs 

which  are  in  phase  with  the  polarisation  wave  induced  by  the  valence 


» 
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electrons. 

However. 

for 

the 

transverse 

component , 

the  plasma 

induced  polarization 

does 

no  t 

add  significantly  to 

the  total 

polarization 

wave. 

This 

is 

because  the 

transverse 

component 

consists  of  the  product  of  a  symmetric  transverse  electric  field 
component  and  an  asymmetric  longitudinal  magnetic  field  component. 
So  its  distribution  is  always  asymmetric  which  yields  only  the 
higher  order  effects.  The  plasma-induced  longitudinal  component, 
however,  is  always  symmetrically  distributed  about  the  center  plane. 
Thus  it  plays  an  important  part  in  the  polarization  wave. 

3«3  Surface  Wave  Generation 

In  this  section,  we  determine  the  intensity  and  distribution  of 
the  EM  field  generated  by  a  nonlinear  polarization  wave  of  phase 
factor  e^^ujt~^kzz^ .  The  most  commonly  used  method  is  a  plane  wave 
approximation  which  assumes  that  both  the  input  and  output  waves  are 
uniform  plane  waves  of  infinite  extent.  However,  this  assumption  is 
only  applicable  to  waves  propagating  within  a  bulk  crystal  of  large 
transverse  dimensions.  For  the  multilayered  planar  structure  in 
which  we  are  interested  the  boundary  conditions  must  be  taken  into 
consideration  and  a  more  rigorous  treatment  is  required. 

The  Green’s  function  method  can  be  used  to  calculate  the 
generated  field  when  boundaries  are  present.  3onsall  and 
Maradudin^  have  reported  the  use  of  a  Green’s  function  formalism  in 
the  presence  of  a  single-interface  structure.  Expressions  derived 
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for  the  generated  field  contain  a  ‘resonant'  denominator  which 
vanishes  as  the  phase-matching  condition  is  achieved.  This  implies 
that  the  generated  field  is  'resonantly  enhanced’  by  the  phase- 
matching  condition.  However,  these  results  can  only  be  obtained  for 
a  loss  free  structure  and  one  which  also  provides  an  unlimited 
propagation  length.  The  manifestation  of  the  'resonance  enhancement 
effect'  arises  from  the  constructive  interference  along  an 
'infinitely  long'  interaction  path.  If  either  loss  or  finite  length 
occurs  then  the  results  have  to  be  modified. 

In  this  work,  a  more  realistic  physical  model  is  used.  The 
structure  is  shown  in  Fig.  12a.  It  is  symmetric  about  the  x*0  plane 
and  has  a  finite  length  L  in  the  direction  of  wave  propagation. 
Absorption  effects  may  also  be  present. 
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The  problem  of  difference  frequency  generation  can  be 
formulated  as  follows.  At  the  difference  frequency  ,  the 

\JT 

macroscopic  field  generated  via  the  induced  polarization  _P  is 

described  by  Maxwell’s  equations: 


VXE  =  -jAwy^H 

(3.3.1  ) 

VXH  =  j  AcjD 

(3.3.2) 

where 

£(r;Au)  =  e0c(x;Acj)n(r;A^)+PNL(r;Au))° 

(3.3.3) 

4 
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r  e(Aw) 


e lx; Au) 


C1 

e2 (Au) 
^  e3 


|x|<t 

t^Mstj+t 

tL+t< |x|<t1+t2+t 


(3.3.4) 


^tl+t2+t 


^  (r  ;  Au) 


'  N  L  -  .  j  (  Awt  -  Ak  z  ) 
P_  (x)eJ  z  J 

0 


|x  [it 

|x  |  >t 


(3-3.5) 


If  the  thickness  of  the  center  layer  is  small  compared  to  the 
wavelength  of  the  polarization  wave  and  all  the  other  transverse 
dimensions  involved,  the  distribution  of  the  nonlinear  polarization 
can  be  approximated  by  a  delta  function  centered  at  the  x*=0  plane. 
That  is, 

PNL(x)  -  PNL  (x)  (3.3.6) 


The  center  layer  can  then  be  eliminated  from  the  analysis.  The 
remaining  problem  is  then  to  determine  the  amplitude  and 
distribution  of  the  field  generated  by  a  nonlinear  polarization 
gource  for  the  simplified  planar  structure  given  in  Fig.  12b.  Two 
different  methods  can  be  used  to  achieve  this  end. 
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3.3.1  Green' b  Function  Foraaliso 

The  Green's  function  formalism  is  widely  used  in  determining 
the  fields  radiated  by  a  given  system  of  sources.  The  Green's 
function  is  the  solution  to  a  given  differential  equation,  together 
with  specified  boundary  conditions,  with  the  source  functions  being 
of  unit  strength  and  localized  at  a  point  in  space.  Then  the  total 
radiated  field  can  be  obtained  by  multiplying  the  Green's  function 
with  the  source  strength  and  integrating  over  the  source 
distribution.  In  other  words,  it  converts  a  differential  equation 
problem  into  an  integration  problem.  In  our  case,  we  have,  by 
combining  eqs.  (  3-3-1 ) —  (  3.3*3),  a  differential  equation: 


.N’L, 


/3xa3xB-6aBV“-6aBko'£Cx;Au))]ES(I;Aa))  =  (k0  /E0)Pa 

(3-3.7) 


where  lc^Aw/c  and  Eg(£,£u)  is  a  component  of  the  macroscopic 
electric  field.  This  inhomogeneous  partial  differential  equation 
can  be  converted  into  an  integral  equation: 

Ea(r;AuO  =  Ea(0)  (r;Aw)-IY(k02/e0)/GaY(Aw|r,r ' '^(r^dV 

(3.3.8) 

where  G  (  Aw  |r,_r’ )  is  a  Green's  function  which  satisfies  the 

I 

equation: 

Z3^3  /3xa3xB'6aBv2‘6aBkoC  (x  ’ Au)^  CBy  (Au),I>l’ )  =  -  <5  (r-r  ' )  6^ 

(3.3.9) 
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and  E^^(_t;Auj)  is  the  solution  for  the  corresponding  homogeneous 


equation. 


To  proceed  further,  it  is  convenient  to  introduce  the  Fourier 
representation  of  (Aid  [r,£‘  )  and  6  (£-£*  ) : 


Gag(  Aid  |r  ,  r  '  )  =  /d^e"-1  z 


-jk  (z-z’) 


gaS(‘k2AwlX*X,)/2TI  (3-3.10) 


6(r,r')  =  6(x-x*)/dk2e'jk2(z'“  '  }  /  2  tt 


(3.3.11) 


The  form  of  this  representation  reflects  the  fact  that  the  system 
characterized  by  the  dielectric  constant  c(xtAu)  (  eq.  (3-3.4))  is 
invariant  in  the  y-z  plane  only.  Substitution  of  eqs.  (3-3-10), 
(3-3.11)  into  eq.  (3-3-9)  yields  the  following  set  of  differential 
equations  for  the  Fourier  coefficient  g  (k„Aid  |x,x' ) : 


e(x;Aid)k02-kz2 


j  kz  3/  3x 


e(x;A<d)ko2+32/3x2-kz2 


jkz3/3X 


U  gf^X-x')^ 

e(x;Aui)k02+32/3x2  (3.3.12) 


The  solution  has  to  satisfy  the  appropriate  boundary  conditions 


it  the  interfaces  as  well  as  at  x=*on. 


I 
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The  Fourier  coefficient  g  of  the  Green's  functions  can  be 

a  8 

obtained  following  reference  (24).  The  unit  source  is  considered  to 
be  localized  at  the  x=0  plane.  The  only  nonzero  Fourier 

coefficients  are  gxx,  Syy»  fizz*  gxz  and  ®zx  which  satisfy  a  set  of 
differential  equations  given  by  eq.  (  3.3.12).  Using  the  boundary 
conditions  given  by  eqs.  A22(a)-(f)  in  reference  (24),  we  solve  for 
these  primary  functions  in  a  straight  forward  manner.  The  results 
are  given  in  Appendix  1.2. 


Now  let's  examine  these  functions  more  carefully.  Functions 
gx„  and  gzz  have  the  same  denominator 


DN1 

»[l+(e2k./e3k2)tank:t2)+(e1k2/e,k1)tank1t1[(c?k  /E-k,)-tank?t,] 
=0  (3.3.13) 


where 


M  ,2 


.  u  2_  v.  2 

K0  e1 ,2 Kz 
-tc02e3  +  kz2 


Both  functions  become  bounded  when  the  denominator  vanishes. 
Similarly,  function  g„x  and  gxx  become  infinite  when 

DN2 

=  (l  +  (£,k./e_k2)  tank^t^]  tank^j*  (c1k^/e,k1)  [tank_,t-,  -  (e^k./e.k^)  ] 

=  °  (3.3.14) 
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infinite  when 


DN3 

= [l+(k./k2)tank2t2]tank1t1+(k2/k1) (tank_,t2-k_/k2) 


Comparing  this  with  the  dispersion  equations  given  in  Appendix 

1.1,  we  find  that  eqs.  (  3.3-0)  and  (  3.3.14)  correspond, 

respectively,  to  the  dispersion  equations  of  the  odd  and  even  TK 

surface  wave  modes  of  the  planar  structure  under  investigation. 

Similarly,  the  root  of  eq.  (  3-3*15)  leads  to  the  even-TE  surface 

wave  modes.  Thus,  the  Fourier  coefficient  g  (k,.^  |x,0)  becomes 

am 

singular  when  kz  approaches  the  propagation  constant  kn  of  some 
surface  wave  mode. 


Combining  eqs.  (  3-3.8)  and  (  3-3-10),  we  can  now  evaluate  the 

components  of  the  generated  field.  Taking  advantage  of  the  delta 

NT 

function  distribution  for  P  ,  the  integral  can  be  reduced  to: 


If  the  structure  is  lossless  and  infinitely  long  in  the  direction  of 


wave  propagation,  as  was  obviously  the  case  in  ref.  (23),  then 
making  use  of  the  fact  that 


/dz'ej  (Vakz)2'/27i  =  6(k_-AkJ 


(3.3.17) 


we  can  carry  out  the  integration  quite  easily.  Each  term  of  the 
summed  series  retains  the  sane  functional  fora  as  the  corresponding 
Fourier  coefficient  g^g(  kzAa)  |x  ,0) ,  except  that  the  integration 
variable  kz  is  now  replaced  by  the  propagation  constant  a  k„  of  the 
polarization  wave.  Consequently  the  amplitude  of  the  generated 


field  demonstrates  the  same  singular  behavior  as  g  ,  i.e.,  it 

aS 

becomes  infinitely  large  when  Ak„  approaches  the  propagation 


constant  kn  of  the  proper  surface  wave  node.  This  is  the  so-called 


phase-matching  condition. 


The  above  conclusion  is  true  only  if  the  interaction  length  is 
infinite  and  no  damping  effect  is  present.  The  violation  of  either 
condition  will  void  the  validity  of  eq.  (  3-3.17)  and  hence  all  the 
results  that  follow.  Unfortunately,  this  is  exactly  what  we  can 
expect  in  the  physical  world.  Therefore,  the  problem  of  EM  wave 
generation  in  the  presence  of  a  planar  waveguide  structure  has  to  be 
re-examined  with  the  finite  interaction  length  and/or  losses  being 
taken  into  consideration. 


Again,  eq.  (  3-3.16)  can  be  used  to  evaluate  the  electric 

fielo  generated  by  the  nonlinear  polarization  source.  However, 

a 
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problem 

arises  with  the 

integration.  Under 

the  condition 

of  a 

finite 

interaction  length 

L,  the  integration 

over  dz'  no 

longer 

yields 

a  delta  function. 

Instead,  we  obtain  a 

sine  function 

which 

has  an  argument  (kz-Akz)L/2.  This  function  describes  a  spatial 

interference  pattern  which  acquires  a  maximum  value  when  kz  is  made 

equal  to  Akz.  In  addition,  if  loss  is  present,  all  wave  vectors 

become  complex  numbers  and  the  expression  for  the  interference 

1  7 

pattern  has  to  be  modified  accordingly 

The  next  step  then  is  to  integrate  over  dk„.  This  takes  into 
account  the  interference  effect  due  to  the  source  distribution  in 
momentum  space.  The  integrand  now  includes  the  Fourier  coefficient 
g  2»  the  sine  function  and  some  phase  factors.  Extra  care  must  be 
taken  in  performing  the  integration  because  of  the  presence  of 
singularity  in  the  function  g  r.  The  integral,  as  discussed  in 

reference  (25),  can  be  decomposed  into  two  parts,  a  branch-cut 
integral  ana  the  residues  contributed  by  the  enclosed  poles.  The 
branch-cut  integral  represents  the  radiation  field  with  a  continuous 
eigenvalue  spectrum  while  the  residues  at  the  enclosed  poles  are  the 
surface  wave  modes. 

It  has  been  noted  previously  that  the  function  g  ^  has  poles  at 

uS 

some  kT  values  which  correspond  to  the  prorogation  constant  of  some 
surface  wave  modes.  For  the  planar  structure  under  investigation, 
there  exist  four  types  of  surface  wave  modes,  i.e.,  ever,  and  odd  TE 
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and  TM  modes.  Generally,  the  excited  surface  wave  is  a  linear 
superposition  of  all  these  modes.  Now  consider  a  special  example  in 
which  the  polarization  source  has  only  an  x-component.  Then  in  eq . 
(  3.3.16),  only  the  terms  associated  with  gxx  and  g2x  give  nonzero 
results.  This  leads  to  the  even-Th  modal  surface  waves.  The 
electric  field  of  the  order  surface  wave  modes  can  be  obtained 

from  eq.  (3-3.16): 

Em  =  (k0^/en)2:;j(  residue  at  k„=k_) 

a  J  0  2  ^  (3-3-16) 


Using  the  expression  given  in  Appendix  1.2  for  g  ,  we  get  the 

a0 

following  explicit  expressions  for  the  electric  field.  In  the  i*^1 
region,  the  electric  field  of  the  order  surface  wave  mode  has 

components : 


(x) 


Exi=-UA/ki)[(dfi  ■  /dx)|DN2=0  ]/[(dDN2/dkJ|k  =k  ] 


E  m=(Akm) [f . 

Z1  11 


z  m 


z  m 


(3.3.19) 


(3-3-20) 


where 

A 


FKC^2eokUi£icos^i 


i  ni 

cosk  t. 


) 


Note  that  k™  represents  the  transverse  wave  numbers  of  the  m^‘  order 
surface  wave  code.  Both  expressions  in  eqs.  (3-3-19)  and  (3.3.20) 
contain  a  factor  F  which  is  the  product  of  a  function 


53 


L  sinc(  (^kz-kn)L/2)  and  a  phase  factor  ej{Au}t-kmz) .  The  s 
function  is  defined  by  sinc(x)=sin(x)/x  where  x=(Ak„-k  )L/2.  It 
negligible  magnitude  for  all  the  eigenvalues  k^  except  the  one  wh 
is  in  the  vicinity  of  Ak„.  Thus  the  surface  mode  which  meets  t 

<u 

phase-matching  condition  becomes  the  dominant  mode.  And  the  fi 
strength  of  the  generated  wave  is  proportional  to  the  interact 
length  L. 

Eqs.  (  3.3.19)  and  (  3-3-20)  represent  the  field  distribut 
of  an  even-TM  modal  surface  wave  excited  by  the  x-component  01 
nonlinear  polarization  source.  Similarly,  the  y  ana  z-components 
a  nonlinear  polarization  source  excite  only  surface  waves  of  even 
and  odd-TM  nature,  respectively.  When  ail  three  components 
present,  only  the  mode  which  satisfies  the  phase-matching  conait 
is  excited. 

3-3.2  Lorentz  Reciprocity  Theorem  Approach 

We  may  also  adept  a  different  approach  utilizing  the  Lore 
Reciprocity  Theorem  in  calculating  the  generated  fields. 
Lorentz  Reciprocity  Theorem  is  applied  in  a  similar  way  to  t 
illustrated  in  ref.  (25).  However,  the  polarization  source  is 
distributed  along  the  direction  of  propagation  as  given  by 
(3.3.6)  and  generates  one  or  more  surface  wave  moces  and  a  radiat 
field.  Since  the  various  surface  wave  modes  are  orthogonal  to 
another,  as  well  as  to  the  radiation  field,  the  total  radiated  fi 
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rr:i rn: 


may  be  represented  as  follows: 


a  1 


1 


for 

the 

forward  direction 

of  propagation 

E  = 

r  „  c  j (Aut-k  z) 
E  a  E  e  v  m 

m  m— m 

+  i 

(3.3.21  ) 

H  = 

E  a  H  ei  (Awt-k  z) 
m  riHn 

(3.3.22) 

for 

the 

reverse  direction 

of  propagation 

E  = 

E  b  E  e-j (Awt-kmt 
m  m— m 

3 

(5-3.23  ) 

H  = 

r  .  u  -j  (Aujt-k  z 
E  b  H  e  J  K  m 

m  m— m 

3 

(5-3.24) 

where  _E _HK  is  the  radiation  field  ana  H_,  the  mth  surface  wave 

mode  field. 

The  total  field  is  a  solution  of  the  equations 


VXE^  =  -jAiop^ 


MI 

VXH  =  j  Aw ( E  E  +  P  ) 


(5.3.25) 


(3.3.26 ; 


while  the  surface  wave  mode  is  a  solution  of  the  source-free  field 
equations.  let  E  and  H_  be  the  total  radiated  field  and  E_,  H_  the 
mth  surface  wave  mode  field  for  the  planar  structure  uncer 
investigation.  If  we  choose  a  closed  loop  S  on  the  x-z  plane, 
consisting  ol  two  transverse  paths  at  and  ^  and  two  longitudinal 
paths  at  infinity,  then  the  2-dimensional  Lorentz  Reciprocity 
Theorem  gives: 
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(3.3.27) 


£  (E  X  H  -  E  X  H)  *nd£  =  ffj AojPNL-E  dA 
S  —  — m  —in  —  J  —  — m 


Following  the  same  argument  as  was  used  in  ref.  (25),  we 
arrived  at  these  results  for  the  amplitude  of  the  radiated  wave  of 
the  m*'*1  mode: 


in  the 


a 

m 


in  the 


b 

m 


forward  direction  of  propagation, 


-lfAwt-k  z)  .  .  _NL ,  , 
JJE  e  m  dxdz 

— m  — 


fe 

Xh  • zdx 

— m  — m 

/E 
— m 

r — <. 

0 

11 

X 

•  A  dNL  j  (k  -Ak  ):, 

■j  AdP^  e  m  z  dz 

fe 

Xh  *zdx 

— m 

— m 

rse 

direc 

tion  of  propagation. 

/E 
— m 

(x=0) 

.  .  N'L  j  (k  +Ak  )z  , 

•j AmPq  eJ  m  z  dz 

fe  Xh  • zdx 

-fli  — m 


(3.3-28) 


(3-3.29) 


where  e^,  are  transverse  electric  and  magnetic  field  components 
associated  with  E_j  and  H^,  respectively. 


In  order  to  show  that  this  method  yields  the  same  result  as 
that  of  the  Green's  function  method,  we  again  consider  the  special 
case  of  having  only  an  2-component  for  the  nonlinear  polarization 
source.  The  dot-product  in  the  numerator  of  eqs.  (3*3.28)  and 
(3*3.29)  suggests  that  the  excited  surface  wave  must  have  a  nonzero 
Ex  component  at  the  x=0  plane.  This  condition  i3  satisfied  by  the 
even-TK  mode  surface  waves.  For  a  structure  of  finite  length  L,  the 
integration  over  dz  again  yields  a  sine  function  which,  in  turn. 
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yields  the  phase-matching  condition,  i.e.,  kn=Ak„,  for  the  nth  order 
surface  wave  mode.  The  amplitude  of  this  phase-matched  surface  wave 
can  then  be  obtained  by  substituting  the  proper  expressions  of  En 
and  2^  into  the  following  equation: 


jAwP  e  (x=0)L 
ox  nxx  ' 

2  /  ”  e  h  cfx 
4  nx  nyH 


(3-3.30) 


where  PQX  is  the  x-component  of  _P 


ML 


obtained  by  performing  an 
integration  over  dx  on  both  sides  of  eq.  (3*3.6).  The  electric 
field  of  the  excited  surface  wave  is  then  given  by 


E  =  a  E  eJU“t-ik  z) 
n-n 


(3-3.31 ) 


Note  that  for  the  even-TM  mode,  the  electric  field  has  only  x  and  z 
components.  And  eq.  (3-3-31  )  gives  the  same  results  as  eas. 
(3-3-19)  and  (3-3-20).  These  are  the  general  expressions  which  are 
applicable  to  both  the  lossy  and  lossless  conditions. 

3-4  Power  Conversion  Efficiency 

3-4.1  The  Phase-Matching  Condition 

The  importance  of  achieving  phase-matched  operation  has  been 
mentioned  previously.  In  the  difference  frequency  generation  case, 
this  requires  the  propagation  constant  (or  phase  velocity)  of  the 
output  wave  to  be  equul  to  that  of  the  driving  wave  of  the  same 
frequency.  Usually  this  will  not  happen  for  intrinsic  GaA3 
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waveguide  as  it  suffers  from  dispersion,  i.e.,  the  phase  velocities 
differ  for  waves  of  different  frequencies.  The  phase  velocity  of 
the  driving  wave  is  in  general  smaller  than  that  of  the  driven  wave. 
In  Fig.  13  we  show  part  of  the  dispersion  curve  of  the  input  TEq 
mode  optical  waves  from  which  the  dispersion  curve  of  the  driving 
wave  is  constructed.  The  result  is  given  in  Fig.  14.  Also  shown 
are  the  dispersion  curves  for  the  driven  MM  waves.  These  are  the 
normal  modes  of  an  intrinsic  waveguide  such  as  that  shown  in  Fig. 
12b,  except  that  all  the  GaAs  layers  are  intrinsic.  In  this 
frequency  range  of  interest,  only  even  TEq  and  TMq  modal  surface 
waves  can  propagate  without  cutoff.  However,  the  dispersion  curves 
of  these  even-mode  driven  waves  always  stay  to  the  left  of  that  of 
the  driving  wave.  That  is,  its  phase  velocity  is  always  larger  than 
that  of  the  driving  wave.  A  phase-matching  operation  is  possible 
only  if  the  KM  wave  can  be  slowec  down  by  some  means. 

Many  ingeneous  phase  matching  techniques  were  reported  which 
used  different  ways  to  compensate  for  the  dispersion  effect.  Ve  are 
particularly  interested  in  the  magneto-optic  phase-matching 
technique4  .  As  we  have  seen  previously,  at  MM  or  SMM  frequencies, 
the  dielectric  property  of  GaAs  material  can  be  controlled  through 
the  plasma  density.  The  plasma  contribution  to  the  dielectric 
constant  is  of  opposite  sign  to  the  lattice  contribution.  The 
degree  of  compensation  can  be  adjusted  by  varying  either  the  plasma 
concentration  or  the  applied  magnetic  field  strength. 
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The  application  of  a  static  magnetic  field  for  phase-matching 


purpose  will  be  discussed  in  Section  3*5.  Here  we  check  the 
feasibility  of  achieving  phase-matched  operation  through  the  control 
of  plasma  concentration.  It  is  illustrated  by  an  example  given  in 
Fig.  15.  For  a  plasma  level  of  ^2.1X10  ^cm  ,  phase  matching  is 
obtained  for  the  driving  and  driven  waves  at  a  frequency  around  270 
GHz  where  their  dispersion  curves  intersect.  A  tuning  range  of  ^210 

GHz  can  be  obtained  by  varying  the  plasma  density  by  %50£.  Both 

27  28 

optical  and  electronic  means  have  been  reported^  ’  for  the  dynamic 
control  of  the  plasma  concentration. 

3.4.2  Absorption  Effect 

The  major  hinderance  in  obtaining  high  efficiency  in  power 
conversion,  besides  the  phase  mismatch,  is  the  absorption  effect  in 
the  waveguide  structure.  There  are  two  major  mechanisms,  the 
dielectric  loss  and  the  free  carrier  loss,  that  cause  attenuation  in 
the  nonlinear  generation  process.  The  dielectric  loss  is  attributed 
to  the  multiphonon  processes  in  the  GaAs  layers  and  the  intervening 
quartz  layers.  In  the  low  temperature  operation,  this  loss 
mechanism  can  be  neglected  for  both  the  optical  and  MM  waves  when 
compared  with  the  free  carrier  absorption  mechanism. 


ft 


» 


ft 


» 


ft  < 


Free  carrier  absorption  occurs  when  the  semiconducting  material  * 

is  extrinsically  doped.  In  an  n-type  GaAs  medium,  electrons  lose 
energy  when  impeded  by  impurity  ions.  Thi3  scattering  process  is 

ft 


ft 
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characterized  by  a  momentum  relaxation  time  or  collision  time  which 

* 

ia  related  to  the  electron  mobility  by  t  "eu/o  .  It  is  well  known 
that  the  electron  mobility  in  semiconductor  material  decreases  as 
the  doping  level  increases.  Therefore,  a  high  impurity 

concentration  will  increase  the  collision  frequency  and  thus  the 
loss.  For  the  input  optical  waves,  this  absorption  effect  is  still 
negligible  if  the  plasma  frequency  is  much  lower  than  the  optical 
wave  frequencies.  This  happens  when  the  plasma  concentration  is 
less  than  10  °cm-  .  For  the  output  KM  waves,  however,  the  presence 
of  an  electron  plasma  brings  an  appreciable  amount  of  attenuation. 
The  collision  frequency  is  comparable  to  the  frequency  of  the 
changing  field  in  spite  of  the  low  doping  level  (  M0  5cn";)  and  low 
ambient  temperature  (  77°K). 

Clearly,  the  consideration  of  phase-matched  operation  and  low 
free  carrier  absorption  presents  two  oppositing  restraints  on  the 
plasma  concentration.  The  phase-matching  operation  demands  the 
slowing  down  of  the  KM  waves  which  is  achieved  by  introducing 
electron  plasmas  into  the  guiding  medium.  However,  since  the 
electron  plasmas  are  accompanied  by  the  presence  of  impurity  ions, 
considerable  loss  arises. 

The  only  way  to  curtail  the  loss  while  at  the  same  time 
maintain  a  phase-matched  operation  is  to  generate  the  electron 
plasmas  by  some  means  other  than  extrinsic  doping.  Thus,  the 


electron  mobility  will  be  maintained  at  a  near-intrinsic  value 
( a/ 2X1  O^cm^/V-sec )  and  the  momentum  relaxation  time  can  be  as  high 
as  20X10  s.  In  Fig.  15,  we  find  the  attenuation  range  from  about 
1-5  to  3*5  dB/mm  for  frequencies  between  160  to  300  GKz.  If 
extrinsic  doping  is  used,  the  electron  mobility  will  retain  only 
half  of  its  intrinsic  value  and  the  attenuation  per  unit  length  will 
be  at  least  doubled. 


3.4.3  Sample  Calculation 

How  we  will  quantatively  estimate  the  efficiency  of  power 
conversion  for  the  proposed  thin  film  structure.  The  efficiency  of 
the  difference  frequency  generation  by  nonlinearly  mixing  two 
optical  waves  is  given  by: 


e  ( Aio ) 


S(Au) 


S(oj1  )  +  S  (  d 2  ) 


(3.4.1  ) 


where  S  stands  for  the  net  energy  flux  per  unit  width  of  the  wave  in 
the  transverse  direction.  At  optical  frequencies,  the  time-averaged 
Poynting  vector  of  a  TE-raode  wave  travelling  in  the  z-direction  is 
given  by: 


<S(x;w)>  =  Re(EXH*)/2 

where  _E  and  H  are  the  electric  and  magnetic  fields. 


(3.4.2) 


The  net  input  energy  flux  S  per  unit  width  is  therefore 
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obtained : 


®(uif2^  a  ,  2  * zcl2C 


t C  —  /“ - ^ - !*<*;«*»,  2)  |2dx]  |l,  2\ 

2f ' 

1  T1  ,2  I1! ,2  I2 


kz (“1  , 2  ^ 


(3.4.3) 


In  the  above  expression  t  gives  the  thickness  of  the  optical 
waveguide  and  e(x;oi^  2)  and  $(x;u)^  2)*  defined  by  eas.  (3.3.4)  and 
(1.1.4),  respectively,  describe  the  dielectric  and  field 
distributions.  The  tern  T1  2  can  be  obtained  from  a  knowledge  of 
the  field  distribution  at  the  respective  frequency.  I1  2  gives  a 
measure  of  the  field  strength.  The  total  power  of  the  input  wave  is 


then  given  by 


=  h  tT  <  oil 


1  ,  2  1  M  ,  2 


where  h  is  the  height  of  the  structure. 


(3.4.4) 


The  next  task  is  to  determine  the  energy  flux  of  the  generated 
output  wave.  At  the  difference  frequency,  the  general  expression 
for  the  energy  flux  per  unit  width  associated  with  an  excited  n**1 
order  surface  wave  node  is  given  by: 

SUw)  -  laj2/”  eriXhn*.Sdx/2  (3.4.5) 

~  OO 

where  _en  and  _hn  represent  the  electric  and  magnetic  field 
distribution  and  an  measures  the  amplitude  of  the  excitation.  Among 


all  the  excited  surface  wave  modes,  only  the  A  branch  of  the  even  TM 


modes  meets  the  phase-matching  requirement.  Therefore,  energy 
transfers  from  the  input  waves  into  the  A  type  MM  wave  mode.  The 
excitation  strength  and  energy  flux  can  be  obtained  by  substituting 
the  proper  expressions  of  the  field  distribution  of  mode  A  into  eqs. 
(3.3.30)  and  (3-4.5),  respectively.  The  results  are: 

-JPoXAu)(Ex1  lx»0^ 

a  -  - - -  (3-4.6) 

2  /” E_ ( x  )  H  ( x  )  d  x 

-00  *  J 

and 

S  (Acj  )  =  |ae  |2/”Ex(x)Hy*(x)dx/2 

=  1  a  |  2  L/ri  E  ■,  H  i*dx^+tc  E  2Hv2*dx+  /"  Hy3*d  x  ] 

0  7  ^  y  \+t2  X3  75  (3.4.7) 

where  (Ex^,  Hxil,  i  =  1 ,3  can  be  obtained  from  eqs. 
(  2.1 .7)— (  2.1.12))  by  replacing  the  thicknesses  T  and  W  with  t1  and 
t2,  respectively.  Note  that  the  strength  of  the  excitation,  ae,  is 
linearly  proportional  to  A10  •  That  is,  the  smaller  the  difference 
frequency,  the  weaker  the  excitation.  It  is  also  linearly 

proportional  to  Pox  which  is  obtained  by  integrating  P^,  the  x- 
component  of  the  nonlinear  polarisation  distribution,  over  dx, 
keeping  in  the  integrand  only  the  contribution  from  the  valence 
electrons.  Thus, 


"e  OX  vkz  1  )  kz  (u  2  ) 

o  " 

1W  2C  0  )c(i;up) 


3  (  *  .  )>3  (  x  ;w  „  ) 


(3-4.8) 
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(3-4-9) 


Pox  "  2/tpxdx  *  tT  1,1  * 

0  v 

Clearly,  PQX  depends  on  the  amplitudes  of  both  input  waves.  And  its 
dependence  on  the  source  distribution  is  summed  up  in  the  terra  T  . 


Now  we  can  calculate  the  energy  flux  per  unit  width  of  the 
output  wave.  By  substituting  eqs.  (  3-4.6)  and  (  3.4.9)  into  eq. 
(  3-4.7),  we  get  the  following  expression: 


SOW)  c  Tp2t2|l1I2l2T3[sinc((A  kz-BA)L/2)]2L2 

(3-4.10) 

where  is  the  propagation  constant  of  the  excited  mode  A  ana  is 
proportional  to  (Acu)  .  Under  phase-matching  operation  and  with  no 
loss,  the  maximum  efficiency  of  power  conversion  is  obtained: 


e(Aw  ) 


tTp2 Ii1i2I2t5l2 

(T,  |I,  \2*T2 \I2  I2) 


(3.4.11) 


It  is  quadratically  proportional  to  the  difference  frequency.  For 
MM  and  SKM  wave  generation,  the  conversion  efficiency  usually  is 
quite  low.  For  example,  if  we  take  L=h=5A=6.3  mm,  two  C02  lasers  of 
output  at  100  V  each  would  generate  a  single-mode  output  wave  at 
A  ^2  50  GHz  with  power  in  the  10-8  W  range.  This  is  of  the  some 
order  of  magnitude  as  the  power  generated  in  a  bulk  crystal  which 
has  a  cross  section  of  1  mmX6-3  mm  and  a  length  of  6.3  mm.  Ideally, 
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a  much  longer  structure  could  be  used  which  will  boost  the 
efficiency  quaaratically  as  the  length  of  the  interaction  path. 
However,  in  practical  cases,  the  free  carrier  loss  has  to  be 


included.  All  wave  numbers  will  then  assume  complex  values  and  the 
sine  function  in  eq.  (3.3*10)  has  to  be  modified  accordingly. 
Similar  to  eq.  (2.29)  of  ref.  (17),  we  also  find  an  exponential  term 
which  severely  shortens  the  effective  interaction  length.  This  is 
the  major  drawback  of  the  use  of  solid  state  plasmas  in  the 
nonlinear  generation  process. 


3-5  Magnetic  Field  Effect 

We  will  devote  the  final  section  to  discuss  qualitatively  the 
role  of  a  d.c.  magnetic  field  in  nonlinear  generation  schemes.  In 
the  presence  of  the  d.c.  magnetic  field,  the  doped  GaAs  medium 
becomes  gyroelectric  and  its  dielectric  property  is  characterised  by 
an  asymmetric  permittivity  tensor.  At  a  moderate  magnetic  field 
strength,  e.g.,  3q=3-8  KG,  the  cyclotron  frequency  is  of  the  order 
of  10  Hz  which  is  much  smaller  than  the  optical  frequency.  Thus, 
the  propagation  characteristic  of  optical  waves  is  not  affected  by 
the  applied  magnetic  field.  However,  for  the  nonlinearly  generated 
difference  frequency  wave,  the  frequency  falls  in  the  same  range  as 
ojc*  Therefore,  cyclotron  resonance  may  play  an  important  role  in 
the  nonlinear  generation  process.  Indeed,  as  was  pointed  out  in 
references  (29)  and  (30),  the  static  magnetic  field  may  have  an 


important  supportive  role  in  the  nonlinear  b 


enavior  of  the 
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conduction  electrons. 


Consider  a  d.c.  magnetic  field  applied  in  the  y-airection 
transverse  to  the  direction  of  wave  propagation  z.  Following  the 
procedure  given  in  Section  3.2.1,  yet  adding  a  d.c.  term,  i.e. ,  Bq, 
to  eq.  (  3-2.6),  we  obtain  the  following  expression  for  the 
nonlinear  polarization  induced  by  two  TEq  mode  optical  waves: 


>NL 


(Aw  ) 


•-O^p 


o  * 

e/m 


OJ  1  O) 


2A(jJ  (  (Aoo )"  '•iiic  ) 


—  I  (Acj'Ococ^X)^(E(w1  )*  E  2  )  )  ]  1 


(3.5.1) 


In  the  above  expression,  the  symbol  X  represents  a  vector  cross- 
proauct  operation.  Clearly,  the  above  equation  shows  a  strong 

dependence  on  the  applied  magnetic  field.  It  exhibits  resonance 
phenomena  as  the  cyclotron  frequency  approaches  the  difference 
frequency.  In  addition,  the  presence  of  the  magnetic  field  causes 
the  coupling  between  the  longitudinal  and  transverse  electric  field 
components  which  may  result  in  an  increase  in  the  nonlinear 

polarizability.  We  will  not  exploit  the  resonance  phenomenon  here 
since  in  this  study  the  magnetic  field  strength  is  predetermined  by 
the  phase-matching  condition  which,  as  will  be  seen  later,  requires 
that  coc-Aio/3.  The  coupling  of  the  longitudinal  ana  transverse 

components,  however,  may  prove  to  be  important  in  the  generation 

process. 
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For  two  input  optical  waves  of  TEg  mode,  eq .  (  3.5-0  yields 
two  nonzero  components.  Comparing  with  eq.  (  3-2.12),  we  find  that 
the  presence  of  a  d.c.  magnetic  field  brings  an  additional  term  to 

MT 

the  transverse  component  of  the  nonlinear  polarization  source  P  . 

cx 

This  term  i3  proportional  to  the  applied  magnetic  field.  It  is  in 
phase  with  the  contribution  from  the  valence  electrons  and,  more 
important,  it  has  a  symmetric  distribution  inside  the  optical 
waveguide.  Therefore,  the  electron  plasmas  now  contribute 
significantly  in  the  nonlinear  generation  process.  Similarly,  the 
presence  of  the  magnetic  field  may  also  introduce  an  extra  term  into 
the  longitucir.al  nonlinear  polarization  component.  However,  the 
induced  polarization  has  an  asymmetric  distribution  inside  the 
waveguide  and  gives  only  a  higher  order  contribution  to  the 
generated  field. 

Now  we  would  like  to  find  the  amplitude  and  distribution  of  the 
MM  wave  generated  by  the  nonlinear  polarization  in  the  presence  of  a 
static  magnetic  field.  We  again  use  the  models  shown  in  Figs.  12a 
and  12b.  The  magnetic  field  can  be  applied  to  the  structure  in  two 
ways.  In  the  first  case,  a  static  magnetic  field  in  the  y-direction 
is  applied  uniformly  across  the  whole  structure.  At  the  difference 
frequency,  the  induced  polarization  nas  two  nonzero  components,  both 
of  which  are  symmetrically  distributed  inside  the  optical  waveguide. 
The  transverse  component  consists  of  contributions  from  both  the 
valence  electrons  and  the  conduction  electrons  while  the 
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longitudinal  component  is  induced  solely  by  the  conducti 
electrons.  In  the  second  case,  the  d.c.  magnetic  field  stays  in  ai 
assumes  opposite  directions  in  the  two  thick  n-GaAs  layers  and  i 
static  magnetic  field  exists  in  the  thin  GaAs  region.  Therefori 
the  conduction  electrons  only  contribute  to  the  longitudin; 

component  of  the  nonlinear  polarisation  while  the  valence  electro! 
are  the  sole  contributor  to  the  transverse  polarisation  componen' 
For  the  generated  MM  wave,  the  presence  of  a  static  magnetic  fiel 
in  the  outer  n-GaAs  layers  changes  these  regions  into  anisotropi 
media.  Thus  the  analysis  given  in  Section  3-3  on  the  surface  wav 
generation  has  to  be  modified. 

We  will  first  study  the  dispersion  characteristic  and  fiel 
distribution  of  the  normal  mode  MM  waves  in  the  presence  c 

magnetoplasmas.  The  next  task  is  to  determine  the  intensity  ar 

mode  structure  for  the  output  waves.  Here  we  do  not  intend  1 

perform  an  exact  analysis  using  dyadic  Green’s  function  or  tl 
modified  Lorentz  Reciprocity  Theorem  on  anisotropic  substrate: 
Rather,  we  will  only  characterize  the  output  wave  qualitatively 
The  output  wave  mode  determination  will  be  based  on  the  followii 
two  criteria.  First,  the  excited  wave  must  have  a  nonzero  electr: 
field  component  along  the  direction  of  the  polarization  sourc: 
Secondly,  the  phase  velocity  of  the  output  wave  must  match  with  th; 
of  the  driving  wave. 
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In  the  presence  of  a  d.c.  magnetic  field,  the  dispersion 
characteristic  of  the  normal  modes  can  be  obtained  from  the 


» 


» 

following  transcendental  equation: 

e-2k2T2  * 

e2k2T2(R1e2k1T1-R2e-2k2T2)-(R3e2k1T1-R4e-2k1Tl) 
e2k2T2(R5e2k1T1-R6e-2k1Tl)-(R7e2k1T1-R8e"2k1Tl)  (3.5.2) 

I 

where 

k  —  - V 

ki  ,3  Y  Ko  '-i  ,3 

p  p  p 

k2  =  -T-k0^e  , 

The  coefficients  R1  to  Rg  are  obtained  from  the  function 

R(  u  ,  k  ,  v  ,  1  ,  m,  n  ,  p  ,  q  )  =  (  k2  +  u  jyn4/C  +ke  gk^/e  5  )  , 

•(k2  +  v  jYn4/5  +  1eeR1  /  e  1  ) 

•(k2  +  mjyn2/r;  +  nE:ek1  /e1  ) 

•(k2+p jYn2/C+qeek5/e5 ) 

in  the  following  way: 
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R6  =  R(-1  1  1-1  1  1  -1  -i  ) 

R?  =  R(  1  -1  -1  -1  1  -1  -1  -1  ) 

Ra  =  R(  1  -1  -1  1  1  1  -1  -1  ) 

where  £  and  n 2  4  are  defined  by  eqs.  (1.2.4)  and  (1.2.6),  with  2,4 
referring  to  each  of  the  two  outer  GaAs  slabs.  If  the  magnetic 

field  is  applied  uniformly  through  the  structure,  we  have  r'i2  =  r'lA 
eq.  (3-5*2).  The  waveguide  modes  can  be  decomposed  into  pure  TE  and 
TM  modes  and  the  propagation  is  reciprocal.  The  field 

distributions,  however,  are  neither  symmetric  or  anti-symmetric 
about  the  center  plane  because  the  structure  is  no  longer  symmetric 
electrically.  The  symmetry  can  be  retained  only  if  the  d.c. 

magnetic  field  assume  opposite  directions  in  the  two  outer  GaAs 
regions.  Then  in  eq.  (3*5.2),  n  2 =  —rl  4  and  the  normal  modes  are 
identical  with  those  obtained  for  the  two  types  of  half  structures 
each  consisting  of  half  of  the  original  structure  and  an  electric  or 
magnetic  wall  at  the  x=0  plane.  If  an  electric  wall  is  in  place, 

the  structure  becomes  an  insulated  image  guide.  Its  TE  and  TM  modes 
have  the  same  dispersison  characteristic  and  field  distribution  as 
the  odd  TE  and  even-TH  modes  of  the  original  structure.  Similarly, 
if  a  magnetic  wall  is  in  place,  the  normal  mode  of  the  half 
structure  correspond  to  the  even  TE  and  odd  TM  modes  of  the  original 
structure. 

According  to  Lorentz  Reciprocity  Theorem,  a  polarization  source 
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will  excite  only  those  waves  of  which  the  electric  field  has  nonzero 
component  along  the  direction  of  the  source  vector.  This  is  still 
true  from  the  modified  Reciprocity  Theorem  on  anisotropic 
structures^1 .  Consequently,  no  TE  node  surface  wave  will  be  excited 
since  its  electric  field  is  perpendicular  to  the  polarization 
source.  The  output  wave  is  thus  pure  TNI  in  nature. 

Now  let's  consider  the  case  in  which  the  d.c.  magnetic  field 
assumes  opposite  directions  in  the  two  thick  GaAs  layers.  The 
dispersion  curve  of  the  odd  TM  modal  surface  waves  has  two  branches, 
S  and  A,  for  both  directions  of  propagation.  However,  surface  waves 
can  only  propagate  at  frequencies  above  some  cutoff  value. 
Throughout  the  frequency  range  of  interest,  waves  of  an  odd  TM 
nature  only  exist  as  evanescent  waves.  The  dispersion 
characteristic  of  the  even  TM  normal  modes  in  the  MM  and  EMM  wave 
range  of  interest  is  given  in  Fig.  16.  A  positive  (negative)  d.c. 
magnetic  field  strength  indicates  the  forward  (reverse)  direction  of 
propagation  and  a  zero  magnetic  field  strength  represents  an 
isotropic  case.  As  was  discussed  in  Section  2.1.1,  there  are, 
generally,  two  dispersion  branches  for  both  directions  of 
propagation.  However,  the  S  branch  wave  is  only  weakly  guided  by 
the  GaAs  region.  Its  dispersion  curve  lies  close  to  the  light  line 
which  is  far  removed  from  the  dispersion  curve  of  the  driving  wave. 
Therefore,  the  phase-match  requirement  cannot  be  satisfied.  The 
dispersion  curve  of  the  A  branch,  however,  moves  away  from  the  light 
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line  toward  high  dielectric  light  lines.  This  is  a  slow  wave  mode. 
Through  proper  control  of  the  strength  of  the  biasing  magnetic 
field,  this  branch  may  be  tuned  to  meet  the  desired  phase-matching 
condition.  From  Fig.  16  we  find  that,  without  the  static  magnetic 
field,  phase  matching  occurs  at  a  frequency  around  270  GHz.  With  a 
moderate  magnetic  field  (  Bq=1 .9  KG  )  applied  parallel  to  either  +y 
or  -$  direction,  a  phase-matched  operation  can  be  maintained  over  a 
bandwidth  of  ~140  GHz. 

Typical  field  distributions  of  the  even-TM  mode  surface  waves 
of  the  canonical  structure  are  illustrated  in  Figs.  17a  and  17b  for 
both  the  forward  and  reverse  directions  of  propagation.  Note  that 
the  electric  field  has  a  tangential  component  which  vanishes  at  the 
center  plane.  Consequently,  there  is  no  coupling  between  the 
tangential  components  of  the  electric  field  and  the  polarization 
source.  On  the  contrary,  there  is  a  strong  coupling  between  the 
normal  components  of  the  electric  field  and  the  polarization  source. 
Therefore,  the  output  wave  is  excited  by  the  normal  component  of  the 
polarization  source  induced  through  the  valence  electrons. 

If  the  d.c.  magnetic  field  is  applied  uniformly  across  the 
structure,  the  TM  modal  waves  no  longer  exhibit  the  same  properties 
as  in  the  previous  case.  As  shown  in  Fig.  18,  the  wave  propagation 
is  reciprocal  and  the  dispersion  curve  of  the  normal  modes  generally 
consists  of  four  branches:  A(,,  Se  and  A0,  SQ.  The  subscripts  e  and 
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o  are  used  to  indicate  the  correspondence  between  these  modes  and 
the  even  and  odb-TM  modes  of  an  electrically  symmetric  dielectric 
waveguide  structure.  Figure  19  shows  the  typical  field  distribution 
of  the  Ag  mode  of  the  canonical  structure.  The  longitudinal  field 
component  no  longer  vanishes  at  the  center  plane.  Therefore  we  can 
expect  a  coupling  between  this  electric  field  component  and  the 
longitudinal  component  of  the  polarisation  wave  induced  through  the 
electron  plasmas.  This  brings  an  increase  in  the  output  power  of 
the  generated  field  due  to  the  additional  nonlinearity  from  the 
conduction  electrons.  However,  this  gain  in  power  is  at  the  expense 
of  the  bandwidth  over  which  the  phase  matching  operation  is 
maintained.  Comparing  Figs.  16  and  18,  we  find  the  bandwidth  is 
greatly  reduced  as  a  result  of  the  external  field. 

More  research  work  has  to  be  done  to  understand  fully  and  take 
advantage  of  the  important  role  of  solid  state  plasmas  in  the 
nonlinear  generation  process.  The  resonance  phenomenon  exhibited  in 
the  nonlinear  behavior  of  magnetoplasmas  may  prove  to  be  essential 
in  overcoming  the  difficulty  imposed  on  the  nonlinear  generation 
scheme  by  the  free  carrier  absorption  effect. 


4.  Conclusion 


Ve  have  exploited  the  feasibility  of  utilising  solid  state 
plasmas  for  potential  MM  ana  SUM  applications,  in  particular,  the 
non- reciprocal  components  and  frequency  conversion.  In  the  first 
part  of  this  work,  the  propagation  characteristics  of 
magnetoplasmons  on  one- dimensional  insulatea  image  guide  structures 
and  two-dimensional  H-guiae  structures  were  investigated.  Both 
showed  strong  non- reciprocity  which  were  usea  as  a  basis  for  the 
design  of  non- reciprocal  devices  sucn  as  isolators  ar.a  differential 
phase  shifters.  Optimisation  of  device  performance  was  applied  on 
the  insulated  image  guide  structures.  Sample  calculations  showed 
quite  acceptable  device  behavior.  The  isolators  were  estimated  to 
give  an  isolation  ratio  of  10:1  ar.a  an  insertion  loss  of  less  than 
0.5  dE/mm  over  a  banawiath  of  65  GHz  in  the  400  GHz  range.  For 
differential  phase  shifters,  phase  shifts  of  35°/mm  to  65°/mm  while 
the  insertion  loss  ox  less  than  1  aE/mm  could  be  obtainea  over  a 
bandwidth  of  ~90  GHz  centered  at  530  GHz.  These  nighly  premising 
results  are  expected  to  encourage  further  experimental  and 
theoretical  work  in  this  direction  to  validate  the  predictions  based 
on  the  current  model  ana  to  develop  models  for  more  sophisticated 
device  designs. 


In  the  second  part  of  this  work,  we  examined  the  process  of 
difference  frequency  generation  of  MM  and  EMM  waves  in  clunar 
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structures  containing  CaAs  material.  Special  attention  was  paid  to 
the  optical  nonlinearity,  phase-matching  condition  ar.d  absorption 
effect  and  their  dependence  on  the  electron  plasma  density  existing 
in  the  structure.  The  presence  of  electronic  plasmas  in  GaAs  layers 
enabled  the  necessary  phase  matching.  However,  it  might  also 
increase  the  attenuation,  thus  reduce  the  power  conversion 
efficiency,  through  the  free  carrier  absorption  effect.  Another 
limiting  factor  for  the  conversion  efficiency  is  its  quadratic 
dependence  on  the  output  frequency.  Thus  for  MM  and  SHM  waves  which 
have  frequencies  far  below  the  optical  frequencies,  the  conversion 
efficiency  stays  low.  A  biasing  magnetic  i'ieia,  however,  may  be 
used  to  overcome  these  aif f iculties.  It  was  found  that  the  presence 
of  a  a.c.  magnetic  field  caused  a  coupling  between  the  ’ or.gitud inal 
and  transverse  electric  field  components  which,  in  turn,  resulted  in 
an  net  increase  in  the  nonlinear  polarisation  induced  through  the 
presence  of  plasmas.  Furthermore,  magnetized  plasmas  also  exhibited 
resonance  phenomenon  in  its  nonlinear  behavior  which  may  prove  to  be 
essential  in  obtaining  a  favored  output. 

More  research  w^ik  is  needed  for  a  quantitative  evaluation  of 
the  nonlinear  generation  on  anisotropic  substrates  containing  solid 
state  plasmas.  This  work  takes  the  first  step  forward  in  that 
uirec  tion. 
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I.  Appendices 


1.1  Normal  Modes  of  a  Symmetric  Dielectric  Waveguide 

For  the  planar  waveguide  structure  shown  in  Fig.  12b,  there 
exist  four  fundamental  modes:  even-TM,  oad-TM,  even-TE  and  odd-TE 
modes.  Each  node  is  a  solution  of  the  dispersion  equation: 


^2 ( ¥■]  ^an^2 *2 ^  +^a (^2+^1  btan^2 ( I •  1  .1  ) 


where  k. 
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In  the  ith  region  the  field  distribution  can  be  found  as  follows: 

TM-mode 

Exi  =  V(ie){x)/(“E0ei) 
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TE-mode 


(x)/aivio 

(d<^-h^(x)/dx)/(jum  ) 
1  0 


(1.1.3) 


where 


<j>-|(x)  =  J(l+b  tan  k2t2)cos  k^cos  k-jX/cos  k-jt-j  for  even  mode 


4>2(x) 
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^1+b  tan  k2t2)cos  k2t,,sin  k^jx|/sin  k-jt-j  for  odd  mode 
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6  (1.1.4) 


b  =,  ^or  ™  mode 

k3/k2  for  TE  mode 
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1.2  Fourier  Coefficient  of  the  Green's  Functions 


Listed  below  are  the  Fourier  coefficients  of  the  Green  s 
functions  for  a  structure  depicted  in  Fig.  12b.  The  unit 
polarization  source  is  located  at  the  x=0  plane  and  the  structure  is 
uniform  in  the  y-direction.  In  the  i**1  region. 


8zz’>=  kifi(x)/(2k0£i  cosk1t1  cosk,t,  DN1) 

gx^  =  ~(jkc/kP  ^ d  5  5  !  }  /  d  x ) 

a(i)=  (jk_D.Nl/k  DN2)g_(_l) 

°  Z  X 

2xx^=  -  jU_dg^^/dx)/kT 

g^)=  hi(x)/(2k1cosk1t1  cosk,t,  DNS) 


(1.2.1  ) 


where 

DN1  =  (l  +  bT2)  +  Ce1k2/e2k1)T],  (b-T,) 

DN2  =  ( 1 +bT2 )  +  (£1k2/£2k1) (T2-b) 

DNS  =  (1+bt2)t1  ♦  (k2/k1) (T2-B) 
b  =  e2k./e3k2  B  =  k  /k, 

T  =  tank  t.  T-,  =  tank^t. 


•  • 


•  • 


(x)  =  [D.Vlcosk^t  ^cosk^  (  x  -t  -  (1+bT.,)  cosk^x]  cosk^t^/sink^t^ 

= [-DNZcos-^t^cosk^ (  x  -t^)+(l+bT0)sink^  x  Jcosk^t^/cosk^t^ 

(1.2.2) 


•  « 


f9(x)  =  sink2(  x  +bcos-7  (  x  -t1-t?) 


f3(x)  =  e'k3C  x  "W 


(1.2.3) 

(1.2.4) 


(x)= [DN’3cosk^t^ sink^ (  x  -t^ -t?)  +  (1  +  BT.,) cosk^xjcosk^t^/cosk^ 

(1.2.5) 


h2(x)=cosk7(  x  -ti -t,) ~Bsink„ (  x  'tj-t.,) 
h_Cx)=e‘k3(  X  'W 


(r .2.6) 

(1.2.7) 


•  « 
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4.  Guiding  and  Control  of  MM-lVaves  by  Surface  Plasmons  - 
Experimental  Feasibility 

In  this  section  we  review  some  of  the  feasibility  estimates  and 
experimental  attempts  which  have  been  made  to  appraise  the  possibili¬ 
ties  of  using  surface  plasma  waves  in  semiconductors  to  guide  MM- 
waves  in  a  practical  environment.  We  consider  here  the  case  where  no 
external  magnetic  fields  are  present  and  where  the  emphasis  is  on  the 
simplest  structure  formed  by  a  single  planar  dielectric  interface. 

There  are,  however,  serious  difficulties  with  diagnostic  evaluation 
of  these  concepts  because  of  present  problems  of  finding  suitably 
flexibile  MM-wave  sources  with  the  semiconductor  structures  of  interest. 
While  our  laboratory  work  so  far  has  yielded  useful  information  about 
practical  design  considerations  and  material  parameters  for  surface 
plasmon  waveguides,  we  have  not  yet  been  successful  with  a  definitive 
demonstration  of  an  unambiguous  propagation  experiment.  Nevertheless, 
with  continuing  advances  in  MM-wave  source  technology,  future  experi¬ 
ments  should  be  possible  with  increasing  versatility. 

4.1  Semiconductor  Material  Considerations 

A  primary  materials  consideration  in  a  practical  surface  plasmon 
MM-wave  device  concerns  the  need  for  low  insertion  losses.  This  implies 
a  requirement  for  a  low  loss  semiconductor  and  a  suitable  geometry  for 
efficient  coupling  (excitation)  of  the  pi  isma  wave  in  a  given  struc¬ 
ture. 

The  attenuation  losses  at  MM-wave  plasma  frequencies  in  an  elemen¬ 
tal  or  a  compound  semiconductor  are  most  likely  to  be  dominated  !•■;.■ 
free  carrier  absorption,  i.c.  from  a  finite'  electron  or  hole 
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With  recent  advances  in  epitaxial  materials  growth  there  are  now 
several  attractive  candidates  to  consider.  Within  the  simple  Drude 
model  for  the  dielectric  response  function  of  an  electron  gas,  we  are 
mainly  interested  in  finding  and  maintaining  a  long  carrier  scattering 
time.  At  the  same  time,  a  small  effective  mass  is  important  in  giving 
the  large  dispersive  contribution  necessary  for  a  reasonable  plasmon 
frequency  at  a  moderate  carrier  density.  In  compound  semiconductors 
such  as  GaAs ,  InSb,  CdTe,  or  HgCdTe,  the  optical  mode  lattice  vibra¬ 
tions  have  a  polar  component  and  this  sets  a  practicable  temperature 
limit  to  a  surface  plasmon  device.  This  follows  from  the  dominance  of 
electron  scattering  by  the  polar  mode  scattering,  typically  at  tempera¬ 
tures  exceeding  77K  in  a  material  of  low  impurity  concentration  (for 

example,  in  InSb,  the  polar  mode  limited  mobility  at  77K  is  approxi- 
4 

mately  10  ) . 

In  the  absence  of  a  significant  phonon  component,  impurity  scat- 
terming  will  play  the  limiting  practical  role  in  the  development  of  low 
loss  surface  plasmon  guided  wave  structures.  Here  we  are  generally 
forced  to  make  an  optimum  compromise  between  an  impurity  concentration 
high  enough  to  yield  the  necessary  carrier  density  but  yet  low  enough 
to  maintain  the  highest  possible  mobility.  For  example,  the  theoreti¬ 
cally  calculated  impurity  scattering  limited  maximum  mobility  in  GaAs 

containing  lO1^  singly  ionized  impurities  per  cm'3  at  77K  is  slightly 
6  7 

less  than  10  cm“/Vsec  (based  on  a  Brooks-ilerring  formulation) .  This 
implies  a  mean  scattering  time  of  4  *  10  ^  sec  for  a  carrier  gas  at 
the  lattice  temperature.  The  precise  free  carrier  absorption  losses 
will,  of  course,  depend  on  the  ratio  of  the  plasmon  frequency  of  the 
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carrier  gas  to  that  of  the  MM-wave  (through  the  spatial  extent  to  which 
the  electric  field  of  the  MM-wave  extends  into  the  semiconductor)  (3) . 

A  rough  estimate  for  the  attenuation  of  the  surface  wave  intensity  for 
a  simple  mode  along  the  interface  can  be  made  through 


L_1  =  2Im(ot)  =  -y" 

A 


3/2 


(4.1) 


where  and  e”  are  the  real  and  imaginary  parts  of  the  dielectric  con¬ 
stant.  The  (free  carrier)  absorption  coefficient  is  related  to 
through  a  =  4IIe'J/2nX,  as  usual. 

Thus  we  may  define  a  crude  figure  of  merit  for  the  case  of  a  single 
dielectric  interface  by  comparing  the  ratios  of  real  and  imaginary  parts 
between  different  semiconductors.  With  the  advances  made  recently  in 
high  quality  growth  of  epitaxial  GaAs,  the  theoretical  limits  of 
impurity  scattering  have  been  nearly  achieved.  Moreover,  there  have 
been  significant  developments  in  the  area  of  MM-wave  integrated  devices 
based  on  this  semiconductor.  Thus,  while  the  calculated  figure  of 
merit  e.g.  for  InSb  is  somewhat  higher,  it  would  appear  that  GaAs  con¬ 
tinues  to  be  the  primary  material  in  the  near  future  for  planned  sur¬ 
face  plasma  wave  guides. 

Recently  it  has  been  shown  that  ultrahigh  electron  mobilities  in 
GaAs  can  be  reached  in  specific  thin  film  structures  of  GaAs,  grown  by 
molecular  beam  epitaxy  (32).  The  idea  here  is  to  utilize  properties 
of  semiconductor  heterojunctions  to  facilitate  spatial  separation  of 
electrons  from  their  parent  impurities.  In  this  way  the  participation 
of  impurity  scattering  has  been  significantly  reduced  in  the  so-called 
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inversion  layers  in  GaAs/Ga^  ^  As  junctions.  Consequently,  low 

x  g  2 

temperature  mobilities  up  to  2  x  10  cm  /Vsec  have  been  achieved 

within  a  high  density  electron  gas  of  bulk  equivalent  density  exceeding 
19  -5 

10  cm  .  For  purposes  of  MM-wave  surface  plasma  propagation  such 
structures  need  to  be  modified  to  allow  a  larger  spatial  extent  of 
the  impurity  ion-free  region.  Typically,  at  (the  rather  long)  MM-wave- 
lengths  one  is  looking  for  plasma  dimensions  transverse  to  propagating 
direction  on  the  order  of  10  micrometers.  Therefore,  the  possible 
exploitation  of  modulation  doping  techniques  would  require  a  multi¬ 
layer  structure  where  a  larger  effective  plasma  volume  would  be  gener¬ 
ated  by  a  multi-quantum  well  composite  of,  say,  several  hundred  layers 
thick.  A  somewhat  different  approach  toward  high  mobility  structures 
has  been  realized  in  the  so-called  nipi  structure  in  which  space  charge 
fields  between  n-i-p  type  junctions  cause  periodic  spatial  modulation 
of  the  conduction  and  valence  band  extrema  (33) .  This  leads  to  electron 
and  hole  separation  in  local  potential  minima.  In  connection  with  the 
MM-wave  surface  plasma  waveguide,  the  nipi  structure  may  offer  another 
attractive  alternative  in  a  multilayer  configuration  which  might  be 
integrable  directly  with  future  GaAs  MM-wave  sources  and  detectors. 

4.2  Excitation  from  an  External  Source 

We  have  made  some  initial  attempts  to  examine  surface  plasma  wave 
characteristics  in  the  laboratory  by  using  far  infrared  lasers  for 
their  excitation.  This  approach  was  chosen  because  of  the  unavailabil¬ 
ity  (to  us)  of  suitable  solid  state  sources.  We  have  encountered  a 
number  of  difficulties  in  developing  an  unambiguous  experiment  in  which 
verification  of  the  propagation  of  an  intrinsic  surface  wave  could  be 
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made.  The  primary  problem  has  been  to  distinguish  between  a  true  sur¬ 
face  wave  and  radiation  which  is  simply  scattered  along  the  surface  in 
the  forward  direction  from  the  input  coupling  optics.  These  problems, 
which  are  discussed  below,  are  not  fundamentally  prohibitive,  but  can 
be  solved  with  an  appropriately  sophisticated  laboratory  approach. 

The  main  source  of  MM-wave  radiation  in  our  work  has  been  a  con¬ 
tinuous  wave  far-infrared  laser,  operating  at  the  (SMM)  wavelength  of 
496  micrometers.  The  active  medium  is  CH^3  gas  whose  vibrational- 
rotational  degrees  of  freedom  are  excited  by  a  pump  CO^  laser.  The  FIR 
laser  was  constructed  according  to  conventional  design  and  consists  of 
a  length  stabilized  resonator  with  a  hole  output  coupler.  Electronic 
feedback  was  employed  in  the  system  to  lock  the  emission  of  the  CO^ 
laser  to  a  position  of  maximum  absorption  in  CH.,  within  the  gain  profile 
of  the  appropriate  line  of  the  laser.  Similarly,  active  feedback  [based 
on  acousto-optic  detection)  controlled  the  resonator  length  and  output 
power  of  the  FIR  laser.  The  spatial  characteristics  of  the  FIR  laser 
were  determined  in  part  by  the  dielectric  waveguide  properties  of  the 
gas  containing  quartz  tube  and  were  of  a  mixture  of  EH-modes.  As  a 
consequence,  the  spatial  characteristics  of  the  output  were  generally 
not  uniform  in  the  transverse  direction.  This  led  to  difficulties  in 
coupling  the  laser  output  to  the  surface  plasma  waves  in  the  semiconduc¬ 
tor  structure,  typically  situated  within  10  to  30  cm  outside  the  laser. 

To  excite  and  detect  the  surface  plasma  wave,  we  set  up  a  prism 
coupling  scheme  similar  to  that  used  earlier  at  much  higher  frequencies 
to  demonstrate  the  characteristics  of  surface  plasmons  in  metals  at 
visible  wavelengths.  (Another  common  alternative,  the  'grating 
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coupling  technique',  was  not  attempted  by  us  here.)  The  basic  idea  is 
to  provide  access  to  the  semiconductor  plasma  which  ordinarily  would 
efficiently  reflect  radiation  at  frequencies  below  the  bulk  plasma  fre¬ 
quency.  A  schematic  arrangement  is  shown  in  Figure  14  where  a  prism 
of  a  sufficiently  high  index  of  refraction  positioned  above  the  semi¬ 
conductor  surface  has  the  role  of  providing  the  appropriately  reduced 
phase  velocity  in  the  direction  parallel  to  the  interface.  In  this 
scheme  ('ATR'  for  attenuated  total  internal  reflection)  a  small  airgap 
is  used  as. a  spacer  layer  between  the  prism  and  the  semiconductor  sur¬ 
face.  The  coupling  from  the  prism  to  the  semiconductor  plasma  medium 
takes  place  across  this  gap  whose  thickness  is  usually  a  fraction  of 
the  wavelength  in  question.  For  a  fixed  incident  frequency  one  now 
varies  the  angle  of  incidence  and  looks  for  a  well-defined  minimum  in 
the  reflectance  from  the  prism  structure  as  a  signature  of  plasma  wave 
excitation.  Analytically,  we  can  derive  the  following  expression  for 
the  reflection  coefficient  for  TH  polarization  in  the  ATR  regime  (the 
subscripts  i,  a,  and  p  refer  to  the  semiconductor  plasma  medium,  the 
airgap,  and  the  prism,  respectively): 


_  1  +  2exp(2ikz.d)Re[r*G(eak:.  -  e.k.^/U^.  *  e.k^)] 
"  1  +  2exp(2ikzid)Re[rpG(eak,.  -  e.k^/C^k,.  ♦  e.k,a)] 


(4.2) 


where  d  is  the  size  of  the  airgap  and  the  complex  quantity  r  is 

rG 

defined  as 
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and  the  wavevector  quantity 
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(4.4) 


where  for  the  incident  angle  a  at  the  prism  base  k  =  n^w/ (c  sin  a). 
The  wavevector  components  and  k^  are  obtained  from  the  disper¬ 
sion  relations 


,  x  ,  2.  2  ,2,1/2  ,  A  ,  2,  2  ,2,1/2  ,. 

k  .  =  ±  (e  to  /c  -  k  )  ;  k  =  ±  (e  to  /c  -  k  )  .  (4.5) 

Z1  1  L  ci  d 


This  set  of  equations  gives  the  necessary  input  for  the  choice  of 
experimental  geometry  and  we  have  used  a  microcomputer  to  calculate 
the  solutions  for  a  number  of  different  materials  as  the  semiconductor 
and  the  prism. 

In  conjunction  with  the  CH.p  496  micron  source,  we  constructed  an 
arrangement  in  which  n-InSb  was  used  as  the  plasma  wave  medium  and 
germanium  as  the  prism  material.  High  quality  bulk  InSb  is  commerci¬ 
ally  available  with  a  range  of  doping  densities  in  the  range  of 
interest.  For  our  initial  purposes,  InSb  was  considered  an  acceptable 
substitute  for  the  generally  more  desirable  GaAs  because  of  cost, 
laboratory  cutting  and  polishing,  etc.  The  choice  of  Ge  was  made 
because  of  its  high  index  of  refraction  and  availability  as  high 
resitivity  material.  Two  equilateral  prisms  were  fabricated  to  facil¬ 
itate  input  and  output  coupling,  respectively.  The  detection  of  FIR 
radiation  was  made  with  a  GaAs  Schottky  barrier  diode  or  a  pyro¬ 
electric  bolometer  and  standard  lock-in  techniques.  The  airgap  between 
the  prism  and  the  polished  InSb  surface  was  chosen  to  be  in  the  range 
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of  50  to  100  microns.  To  verify  the  action  of  the  prism  as  an  input 
coupler,  the  reflection  coefficient  was  measured  for  each  prism  for 
InSb  at  room  temperature.  A  distinct  dip  was  seen  in  the  reflectance 
R(a)  over  a  narrow  range  of  values  for  the  external  angle  of  incidence 
a  with  respect  to  the  prism  side.  The  amplitude  of  the  dip  was 
usually  rather  small,  however. 

In  order  to  study  the  propagation  characteristics  of  a  surface 
plasma  wave  in  InSb,  the  two  3  mm  wide  prisms  were  mounted  with  spacers 
on  the  flat  surface  of  a  polished  wafer  of  InSb,  with  a  variable 
separation  between  them  typically  adjusted  to  about  1  cm.  The  entire 
arrangement  was  then  mounted  into  the  coldfinger  of  a  liquid  nitrogen 
cooled  dewar  at  a  nominal  temperature.  (The  cooling  is  necessary  to 
reduce  the  attenuation  by  phonon  assisted  free  carrier  absorption,  as 
discussed  above.)  The  FIR  laser  radiation  was  focussed  onto  the  input 
coupling  prism  with  a  TPF  lens.  From  this  point  on,  however,  serious 
difficulties  were  encountered  in  attempts  to  make  unambiguous  observa¬ 
tions  of  the  plasma  wave  propagation.  Invariably,  signals  with  a 
reasonable  noise  ration  were  detected  off  the  output  coupler  but  it 
was  unclear  whether  they  were  originating  from  true  surface  modes  or 
whether  they  represented  FIR  radiation  which  was  simply  diffracted 
into  grazing  along  the  surface.  The  origin  of  these  ambiguities  was 
determined  to  be  our  laser  source  whose  spatial  beam  profile  was 
quite  nonoptimal  for  these  purposes.  (The  problem  of  spatial  non¬ 
uniformity  of  waveguide  FIR  lasers  is  a  well-known  drawback  of  these 
otherwise  relatively  simple  devices.)  A  good  deal  of  effort  was 
expended  in  attempts  to  influence  selectively  what  should  have  been 


» 

a  proper  surface  mode  signal,  but  with  little  success.  For  example,  ~ 

placing  a  beam  block  on  the  InSb  surface  strongly  attenuated  the 

output  signals.  Since  a  sizable  fraction  of  the  surface  plasma  wave 

field  was  calculated  to  extend  into  air,  however,  such  a  check  did  - 

not  yield  a  clear  result. 

From  these  initial  experimental  attempts  it  is  clear  that  much 
care  needs  to  be  exercised  in  the  design  of  a  proper  diagnostic  anu  ^ 

evaluative  experimental  program  to  demonstrate  the  feasibility  of 
practical  surface  plasma  MM  or  SMM-wave  guiding.  Much  consideration 
needs  to  be  given  to  a  versatile  (preferably  frequency  tunable)  - 

spatially  well  collimated  and  coherent  source  whose  output  can  be 
directed  into  the  chosen  structure  with  small  angular  uncertainty. 

The  air/semiconductor  interface  which  we  chose  to  examine  may  not  be  ' 

the  best  choice  for  a  waveguide  either  (although  it  is  certainly  the 
simplest).  A  slab  geometry  could  offer  an  advantage  in  permitting  the 
input  and  output  couplers  to  reside  on  the  opposite  sides  of  the  semi-  ^ 

conductor  wafer.  Finally,  with  the  emergence  of  multilayer  epitaxial 
structures  it  may  be  well  worthwhile  to  consider  the  location  of  the 
surface  active  medium  as  a  'buried'  layer.  The  coupling  of  the  desired  ^ 

MM-wave  signals  may  then  be  possible  by  e.g.  grating  techniques  where 
the  corrugation  itself  would  be  buried  at  or  near  the  suitable  inter¬ 
face. 
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Figure  2. 

Effective  dielectric  constant  for  anisotropic  n-GaAs 
with  d.c.  magnetic  field  perpendicular  to  the  propagation 
direction. 
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Figure  3.  Dispersion  curves  for  the  lowest  order  single  air-GaAs 

interface  and  single  dielectric-GaAs  interface.  Long-dashed 
lines  define  radiation  regions  where  <  0. 
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Figure  5a. 

Dispersion  diagram  for  insulated  image  guide  structures 
with  up  3  10'3/s,  uc  =  10l2/s,  t  3  8  x  lcH2s,  c(°)  =  13, 
c,  =  4,  C3  =  1 ,  T  =  25um  and  W  =  50um.  Forward  direction 
Ci  propagation. 
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Figure  6a. 

Isolation  ratio  and  insertion  loss  as  functions  of  frequency 
for  three  GaAs  slab  widths.  The  other  parameters  for  the 
structure  (Figure  1)  are  ej  =  4,  e  =  1  and  T  =  25um. 
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Figure  7.  Differential  phase  shift  and  insertion  loss  as  functions 
of  frequency  for  three  device  configurations.  Here, 
ci  =  13  and  cq  =  1  (see  Fig.  1). 
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Figure  8.  Differential  phase  shift-insertion  loss 
ratio  and  insertion  loss  as  functions 
of  frequency  for  the  sane  cases  as  in 
Figure  7 . 
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Figure  9.  H-guide  structure.  Wave  propagates  along  z-direction 


Figure  10.  * 

Regions  of  'true'  and  'generalized'  surface  waves  defined 
as  functions  of  normalized  frequency  and  wave  number. 
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Figure  lib. 

Dispersion  curves  of  single  interface  and  asymmetrically 
loaded  slab  structure  with  e,  *  4,  c,  *  1,  H  «  0.4(2::c/u  ) 
and  k  c/Wp  =  0.24.  1  J  P 
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Figure  13.  Dispersion  diagram  of  the  input  optical  waves 
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Figure  15.  Dispersion  curves  of  driving  and  driven  waves  in  planar 
intrinsic  GaAs  structure. 


107 


:ield  distribution  of  Even-TM  modes 
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